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ABSTRACT 


It is difficult to propagate a diffraction-limited laser 
beam through the atmosphere, since the atmosphere contains 
random index of refraction fluctuations. Two parameters that 
characterize the atmosphere for optical propagation are the 
atmospheric isoplanatic angle, Be and the refractive 
turbulence structure parameter, cus This dissertation deals 
with improved methods for measuring SPA and с profiles using 
optical techniques. 

By apodizing the receiver telescope aperture, one can 
improve the weighting function for isoplanatic angle 
measurement substantially over previous systems. We find 
that the weighting function is not significantly affected by 
inner scale changes with altitude and that the error in 
isoplanatic angle measurement from strong low altitude 
turbulence (z < 1 km) with this weighting function is small. 
Data collected with the improved isoplanometer shows temporal 
trends in the isoplanatic angle on the order of 90 seconds 
that have not been observed before. 

Direct inversion of the amplitude covariance function 
(including aperture averaging effects) to yield refractive 
turbulence profiles is known to be ill-posed. I suppress 
this condition using Tikhonov regularization and reproduce 
refractive turbulence profiles from actual ə data vith some 


Success. 
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I. INTRODUCTION 


If one wishes to propagate a spatially and temporally 
coherent electromagnetic wave through the atmosphere or any 
random medium without distortion, one must somehow 
compensate for the random phase and amplitude perturbations 
induced by the medium. Adaptive optics is the field of 
physics that attempts to compensate for propagation through 
the atmosphere, principally by deforming mirrors and 
unconventionally by using nonlinear optical materials. In 
order to successfully conjugate the effects of the 
atmosphere, knowledge of certain atmospheric parameters is 
crucial. These parameters are the isoplanatic angle (9,), 
the spatial coherence length of the atmosphere, (r,,; og) anc 
vertical profiles of the refractive turbulence structure 
parameter, PEA This dissertation deals with the theory and 
measurement of the isoplanatic angle and the refractive 
turbulence structure parameter. Currently, reliable systems 
exist for the measurement of the spatial coherence length of 
the atmosphere [Refs. 1, 2]. 

The isoplanatic angle is an angular measure of spatial 
coherence in the atmosphere. Walters of [Ref. 1] has 


developed three generations of instruments (known as 


isoplanometers) to measure the isoplanatic angle. This work 


has to do in part with the development and test of a fourth 
generation isoplanometer having an improved weighting 
function. This instrument is optimized for night use and 
provides temporally high resolution data (1 second samples), 
able for the first time to measure trends with periods less 
than five minutes. Two fourth generation isoplanometers are 
currently operating at different locations. 

Of the three parameters noted, the most important is the 
refractive turbulence que profile. Both the isoplanatic 
angle, D and spatial coherence length, r, are functions of 
ә Hence, knowledge of the вы profile not only gives a 
picture of turbulence with altitude, but also gives the 
measures of spatial coherence, lo and E Unfortunately, 
high vertical resolution profiles of pe are difficult to 
measure. Several measurement techniques exist, each having 
its strengths and weaknesses. These techniques will be 
discussed later. This research attempts remote measurement 
of Ss by direct inversion of the stellar scintillation 
amplitude covariance function. Mathematically, the problem 
is ill-posed meaning that noise in the data makes inversion 
impossible. Different regularization techniques exist to 
make the problem tractable. This work is an extension of 
[Refs. 3, 4] using the Tikhonov regularization technique 
[Ref. 5] to invert the covariance function for two finite 


but arbitrary equal apertures. 
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II. BACKGROUND 


A. GENERAL DESCRIPTION OF TURBULENCE 

Turbulence has certain characteristics that differen- 
tiate it from other flow patterns. In [Ref. 8], Lumley 
briefly describes these characteristics as I will do here. 
The first quality of turbulence almost goes without saying; 
turbulence is irregular. This "quality" has a profound 
effect on how one approaches a mathematical model or 
description of turbulence. Fluid dynamics for years relied 
on the application of Newton's laws to fluids. With the 
Navier-Stokes equations in hand, one in principle should 
be able to solve any problem. However, in the case of 
irregular flows, the randomness makes this approach 
impotent and we are left with having to use statistical 
methods. 

Diffusivity is a property responsible for rapidly 
mixing and spreading basic physical quantities like 
momentum, mass and heat. This property is present in all 
turbulent flows. 

Turbulent flows are always characterized by large 


Reynolds numbers. The Reynolds number, Re, is given by 


Ji 


where v is the velocity, L is a characteristic dimension of 
the entire flow and v is the viscosity. As the Reynolds 
number moves beyond a certain critical value, Rec, 
(nominally Rec ~ 100, for pipe flow) [Ref. 9], the flow 
becomes unstable giving way to turbulence. The value of the 
critical Reynolds number depends on the type of flow. Re, 
for the atmosphere is approximately 2000 [Ref. 14]. 

Another important property of turbulence is that it is 
characterized by three dimensional vorticity fluctuations. 
The vorticity fluctuations find their origin in the velocity 
field. In order to support the random vorticity fluctua- 
tions, the random velocity field must also be three 
dimensional.  Lumley's example of a two dimensional 
atmospheric cyclone is worth noting. Large scale cyclonic 
behavior responsible for weather patterns is itself not 
turbulent. This does not mean, however, that the cyclone is 
not made up of smaller scale turbulent bodies. 

Within every turbulent flow is a mechanism that 
dissipates energy. The dissipation occurs at small scale 
sizes and is a function of the viscosity of the fluid. The 
dissipated energy is internalized into the fluid typically 
as heat. Clearly, if one increases the rate at which energy 
is added to a fluid, the dissipation rate, c, must also rise 
proportionately to conserve energy. This occurs after the 


energy increase has reached the smallest scale sizes. 
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The last two properties in Lumley's list are rather 
obvious but still should be noted. First, turbulence lies 
within the continuum model of fluid mechanics. The scale 
sizes in turbulence are typically much larger than molecular 
sizes. Thus, molecular diffusion does not significantly 
affect turbulent flows. And lastly, turbulent flows in 
fluids are indeed flows. In principle, as I mentioned 
previously, the Navier-Stokes equations provide the general 
equations for fluid flows. For turbulence, one is concerned 
with the statistics of the Navier-Stokes equations. 
Unfortunately, the solution to this set of equations is not 
known since the number of unknowns exceeds the number of 
equations. This is known as the "closure" problem of the 
Navier-Stokes equations  [Ref. 10]. 

As was mentioned earlier, turbulence is an irregular or 
random process. Finding exact scale sizes, dissipation 
rates, etc., is virtually impossible. However, it is 
possible to extract order of magnitude values by dimensional 
analysis. This technique is used in the turbulence field to 
mathematically model the phenomenon. Other techniques 
include asymptotic invariance (the model should behave 
properly as the Reynolds number approaches infinity) and 
local invariance (depending on the time and length scales in 
the turbulence, one may be able to assume that the 


turbulence is everywhere similar). These three techniques, 


|. 


in varying degrees, find their way into the treatment of 
turbulence. 

All mechanical processes must have an energy source. In 
turbulence, the source is principally shear flows within the 
mean flow. This is true for turbulence outside the 
atmospheric boundary layer. The boundary layer is formed 
by the interaction of the flow and the earth's surface given 
the no slip condition [Ref. 8]. Buoyancy can also be an 
energy source. With an energy source available, a laminar 
flow can become turbulent. The laminar portion of the flow 
is modeled by the linearized theories valid for small 
perturbations much like the classical or quantum harmonic 
oscillator. The turbulent portion of the flow is mathe- 
matically treated using an asymptotic theory valid for high 
Reynolds numbers. This inevitably leads to chaos in the 
transition region. 

In turbulence there are a number of different length 
scales. A description of all the different length scales is 
a key component of similarity theory [Ref. 10]. However, 
two important ones are covered now. First, the outer scale, 
typically noted as Loe is the size of the turbulence at 
onset (on the order of meters). The outer scale, given by 
Lo - Re v/v, is a function of the Reynolds number, viscosity 
and velocity of the fluid as mentioned earlier. In the 


outer scale, viscosity does not play a role in the 
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dissipation of energy since inertial terms dominate the 
equations of motion. However, at the inner scale or 
microscale (on the order of millimeters), viscosity 
internalizes the energy of the turbulence as heat. Clearly, 
from what has already been said, the key parameters at the 
inner scale are the dissipation rate and viscosity. Using 
dimensional analysis, Kolmogorov formed length, time and 
velocity scales valid at the microscale [Ref. 11]. These 
are called the Kolmogorov microscales and their exact 
expressions can be noted in [Ref. 8]. 

Now, we can construct a physical picture of the 
formation and dissipation of turbulence. The onset of 
turbulence occurs at the outer scale. If the Reynolds 
number at the outer scale is again greater than the critical 
Reynolds number for turbulence formation, the turbulence at 
the outer scale will again break down. The formation of 
continually smaller eddies occurs until Re < Re, and the 
turbulence is at the inner scale. Hence, energy from the 
ordered flow moves through the turbulence "cascade" to the 
inner scale where it is dissipated in random processes. The 
region of scale sizes between the inner and outer scale is 
known as the inertial subrange. The energy transfer from 
one scale size to the next is adiabatic in the inertial 
subrange. Hence, the entropy increase occurs at the inner 


scale. This view of turbulence is appealing since it is so 


INS 


closely related to the second law of thermodynamics. We 
might expect that the inertial subrange has a well defined 


Spectrum. It does, and this will be covered later. 


B. STATISTICAL DESCRIPTION OF TURBULENCE 
1. Random Variables 

Since turbulence is irregular and random, it is most 
effectively described statistically. Hence, variables like 
the velocity v(t), become random variables having a mean and 
higher statistical moments. In many stochastic processes, 
the statistical moments do not change with time or are 
stationary. However, with atmospheric turbulence this is 
not the case. The nonstationarity of turbulence makes the 
meaning of long term averages, etc., questionable. In 
addition, the replacement of ensemble averages by time 
averages (ergodic assumption) presents additional 
difficulties. This complicates the description of 
turbulence considerably. 

2.  Homogeneity and Isotropy 

Along with stationarity, most well-behaved random 
variables are homogeneous and isotropic.  Homogeneity means 
that statistical quantities do not change with a Galilean 
transformation of coordinates. This implies that 
mathematically, functions describing statistics at Tr, and r, 
depend only on the difference 112 = Ly - T2. Isotropy 


implies a symmetry in rotations of Y} and r5. Hence, if the 
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statistics are homogeneous and isotropic, functions 
depending on [ai and To can be represented by İr) = Fol. 

Unfortunately, the atmosphere provides neither 
homogeneous nor isotropic random variables. However, 
Tatarski in [Ref. 6] adopts the idea of local homogeneity 
and isotropy. Over some region R, comparable to the outer 
scale Lo! the statistics of the random variables do not 
change with translations and rotations. 

In a similar manner, Tatarski handles the non- 
stationary problem. He notes that the functions can be 
stationary over time increments or stationary increments. 
Hence, it is possible to tolerate slow drifts in the mean 
and other moments. 

ori itructure Functions 

Another way to cope with the non-stationary problem 

is to define a function in terms of a difference. Tatarski 


defines the structure function for temperature as 
IE) EA (25 1j 
әз 2 1 
< > denotes an ensemble average. We should note that the 
structure function depends upon the vectors Ii and T3. PE 


we assume homogeneity and isotropy, this reduces to 


D = r -r D) oe OB 
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where r - |r. - ril. Also, in differencing T(r,) and T(r.), 
effects of a changing mean are removed. 

Kolmogorov showed [Refs. 6, 8] by dimensional 
analysis that the structure function for temperature 


fluctuations in the inertial subrange is given by 


22/3 


T QUAS 


Dn (r) = C 


where SCH is the temperature structure parameter. In EM 
propagation, one is more concerned with index of refraction 
than temperature. So, we can define an index of refraction 


structure function, D (r); similar to (2.2) and (2.3), 


2,2/3 


A (2.4) 


D (r) = C 

where Se is the refractive turbulence structure parameter. 
In order to model the turbulence correctly, one must use 

quantities that are not affected by position in the fluid. 
We refer to such quantities as "passive additives". 
Temperature is not a passive additive since it changes with 
vertical displacements. However, potential temperature, 9, 
given by 6 - T - rz where T is the adiabatic lapse rate and 
Z is altitude, is a passive additive since it is corrected 
for changes in altitude. 

Figure 2.1 taken from [Ref. 14] shows the behavior 


O Ë Dj (r) as a function of log(r). Clearly, if r < 1, the 


ES 


5 SX103409A OM] UOOM39Og 
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(14-41) Bol 
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Lu Bol 
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2 
inner scale, Dip (1) = 0 and for r > Lo” DZ (r) approaches 20, 
asymptotically as the fluctuations become uncorrelated. 
Hence, the structure function has a well defined slope in 


the inertial subrange. 


4. Covariance Functions and Spectra 


Along with structure functions, covariance functions 
(or correlation functions if normalized) and spectra are 
useful in describing random processes.  Conveniently, 
covariance functions and spectra form a Fourier transform 
pair. The covariance between two random variables S and T 
is given by 

Bor = SUD nM B «T (r4) 5] [S (r2) = <S(r.)>]>, (22500 
where the average value of T, <T(r)> has different defini- 


tions if T is continuous or discrete. Often one uses the 


autocovariance function 


Buet = «IT (r4) - SEL) ^D GEO - <T(r5)>]>. (2.52) 


If the random variable T is homogeneous or locally 


homogeneous Вр- becomes 


Bop = <T(r,)T(r,)> - SES EN : (2.6) 
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If we choose <T>=0 this reduces to 


Bop = <T(r;)T(r.)> . (2m 7) 
There is a relationship between structure functions and 


covariance functions given by, [Ref. 6], 
Dm(r) = 2[Bmm(0) - Bmm(r)] - (2-8) 


As mentioned earlier, the covariance function and 
the power spectral density form a Fourier transform pair. 
The power spectral density gives one insight into how 
different frequency components contribute to the variance. 


For one dimensional spectra, the transform pair is 


c 


J e Ik B(r) ar (2. 9) 


= o 


W(K) 


and 


oo 


= f DAMIEN. (As LO) 


— OO 


Br) 


But, before these relations hold, we must realize again 


that the functions involved are random. A random function 
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f(r) can be decomposed using Fourier-Stieltjes integrals 


given by 


f(r) = yes do (K) (2-01 


where d$(K) is a random complex amplitude. The Fourier- 
Stieltjes integral must be used since the random 
fluctuations are uncorrelated over non-overlapping intervals 
(Ref. 12]. Hence, f(r) does not have a derivative and is 
non-Riemann integrable. 

With the substitution of the Fourier-Stieltjes 


integral, the covariance function is 
Fi Kiril иси 4 
B(r;tr,r,) = € <A DR) > x; (2.129 


If the medium is homogeneous, B(ri*r, гу) = В(г), since an 
average must not depend on the location within the medium. 
For the double integral to only depend on r, the following 


must be true, 
* 
<do(K)do (K')> = 6&(K-K')W(K) dKdK' S (2.159 


We note that «de (K) de  (K')» is in the form of a mutual 


coherence function (MCF), which represents the overlap of 
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dé (K) and ü gön Since the do's are uncorrelated except at 
K = K', we get the Dirac delta function with the spectrum 
(BORG ving the gamplsitudex= Substituting 42.43) into (2.12) 
and integrating over K' gives (2.10). 

Since we are dealing with fluctuations of the random 


variable T, the mean value is assumed zero giving 


on = f W(K)dK . (2. 14) 


From (2.8) and (2.10) we see immediately that 


со 


fu - cos(rK)] W(K)dK , (2-15) 


= 00 


D(r) = 


alt 


since B(r) is an even function. 
Tatarski [Ref. 6] develops an expression for W(K) 
based upon the form of the structure function for 


turbulence. The one-dimensional spectra W(K) is given by 


W(K) = a 


, 


where a is a constant and e is the dissipation rate. 
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Vector fields like the velocity, v(r), are expressed 
in three dimensions. Hence, their spectra are three 
dimensional. 


We will let (K) denote this power spectral 


density. (K) is given by 


+K) = || are ас | eg Tee 


The covariance function is defined similarly.  B(r) is given 
by 


a —— 5 Joer? e(K) . (2.18) 
(25) 


Equation (2.18) 





is also developed from a Fourier-Stieltjes 


approach requiring homogeneity. For spherical coordinates, 


djr = r” sinedrdoede 


, 


- 


and the @ and 6 integrals in (2.17) can be completed giving 


oo 


Ф(К) = : J B(r) sin(rK)rdar "`". (2.19) 
21 K 0 





A similar relation for the covariance function can also be 


found. Tatarski uses the relation 
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Ф(К) = == (2.20) 


to evaluate the three dimensional spectrum in terms of the 
one dimensional spectrum. This gives 


24711/3 (2.21) 


Ф(К) = aC» 
for the passive additive, index of refraction, where a is a 
constant. This is known as the Kolmogorov spectrum (since 
Kolmogorov did the original work), valid in the inertial 
subrange. Spectra have also been developed for the pro- 


duction and dissipation regions of turbulence [Ref. 13]. 


C. EM PROPAGATION THEORY IN TURBULENCE 
Clifford [Ref. 15] covers Tatarski's work [Ref. 6] for 

weak turbulence. This treatment is analogous to any 

application of first order perturbation theory. The 

following is Clifford's treatment in a more condensed form. 
l. Solving the Wave Equation 

lot 


Consider a sinusoidal wave of the form e where 


the atmosphere has zero conductivity and unit magnetic 
permeability. For these conditions, Maxwell's equation in 


Gaussian units become: 


v ° H =O, (2322) 
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vx E = ikH, (223) 


or (2.24) 


< 
X 
m 

Il 

| 
вә, 
= 
O 
D 


V ° (n”E) - Ü à (29825) 


where k is the wavenumber of the electromagnetic wave and n 
is the index of refraction. We generate the wave equation 
for E in the usual way by substituting (2.24) into the curl 


of (2.23) ugiiwiimg 


-V”E + v(v ə E) = K> . (2.26) 


We can expand (2.25) and solve for V * E. Substituting into 
(2.26) gives 


A A O (ZZ) 


If we assume that depolarization effects are negligible 
(i.e. assume isotropy), the last term of (2.27) can be 
neglected. Clearly, one need only consider a single scalar 
equation, since all the components follow in the 


mathematics. Thus, the scalar wave equation becomes 


V Eocto konmmbo- 1. (2.285 
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Clifford now applies the method of small pertur- 
bations. This means that E can be expressed as a series of 


decreasing values, 
исси ки : (20 


EG corresponds to the unscattered wave, E, to single 
scattering, E, to double scattering and so on. In some 
cons the series may diverge, particularly if ES is 
large or the path length long. This leads to the phenomenon 
of saturation. Clifford uses the first order approximation 


giving 


and 


WEE tok BE Fo ZKPN.E, 509% (2:30) 


All terms of order two have been neglected. Clifford like 
Tatarski assumes a unit amplitude wave propagating in the z 
direction. This allows us to use the standard Green's 
function solution for E, (r), 

ik|r-r'| 


55. DES F: 2: ne... 25 
Eau) = де d'r E [2k n, (r')e Мә 2 31) 


Ҹо1 ə. 
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For laser propagation, most of the monochromatic light will 
be scattered into a small forward cone. Hence, backscatter 
is negligible. Clifford now applies the Fresnel 


approximation to E (T) giving 


2 ikz / š Lu dd ne FA 
— ke B=, ТК | (p 0 D] l 


where p are the coordinates transverse to the propagation 
direction z. 
2. The Rytov Approximation 
At this point both Clifford and Tatarski apply the 
Rytov approximation which assumes that the solution of the 


stochastic wave equation is of the form 

E = Ae ; (2-338 
where A is an amplitude and s is a complex phase. The Rytov 
solution involves the perturbed wave equation (2.30) divided 


by 50: Tatarski gives this equation as 


IE + kn” (6) = log A ss 


This combination of solution and wave equation is known as 


the method of smooth perturbations. The name comes from the 
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fact that we require the perturbations to be small over the 
distance of a wavelength. This condition is less restrict- 


ing than those for the method of small perturbations. 


Brom (2433). the ratio of E/Eg where E = Eg + Es 1s 
E 

Би Ј А | 8 

[m = 1 + En = Ay expfi(S So)] ° (265 


Taking the natural logarithm of (2.35) and expanding in a 
pover series gives Clifford"s result, 


E A 
SÉ : (2. 30) 


Thus, the Rytov approximation has allowed us to obtain 
amplitude and phase information from the Green's function 
solution to the wave equation.  Tatarski and others have 
developed the convention 


I 
x = &n(A/Ag) = — E 


O 


and 


This is convenient since an experiment measures the log 
amplitude x. It is trivial to obtain the forms for Ај/А, 


and (S - So) from (2.32). But it is not so trivial to move 
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on from there. The problem is completing the integral over 
the random function n} (r). If one considers the plane wave 
propagating in a layered media, where at a given z, there 
exists statistical homogeneity, we can expand n} (T) in a 
two-dimensional Fourier-Stieltjes integral. We also assume 
that each layer is uncorrelated so that an integral along 
the path is a sum of incoherent additions. In the notation 


of Clitterd, n4 (r4!) becomes 
nj (r') = IE T? , (2 408) 


where dv (K,z") is the complex Fourier-Stieltjes amplitude. 


Substituting this into the relations for A,/Ag and (S - S.) 


gives 
— Z 
5 - k^ dz' canon (Ky zw) aiser? COS cos a er GI 
S: (r) T (2-2) sin cos a er 277) 
0 
(2. 39) 
where d?r" = а2 "421". x(r) refers to the cosine of the 


bracketed quantity while S} (T) to the sine. One can perform 
the p' integrations leaving the Fourier-Stieltjes integra- 
tion and the z' integration. Now we can construct the 
two-dimensional covariance functions exactly as before. 


B_ (0,2) is 


— m — k oR 
B, (6,z) € cx(9) £ 9, z) x (9),z)? . (2.40) 
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Tatarski finds the spectrum from (2.40) to be 


L 
2 
tio WK = = 2 . WD K (L-z) 
ЕР, (К,2 2"=0) PY (K, 0) 27k Por sin ҹ az 
0 
(2.41) 
and 
: 2 
a 2 2 K (L-2) 
F_(K,0) = 21k Por COS EE dac 
0 


where L is the path length and (K) is the three-dimensional 
Kolmogorov spectrum. 

oe. be AB ygens=Fresnel Approach 

As we have seen, Tatarski and Clifford use the 

differential equations approach to solving the turbulence 
problem. Since we have linear operators, etc., the 
identical problem can be formulated in terms of integral 
equations. This is Huygens-Fresnel theory applied to 
propagation in turbulence. Lutomirski and Yura [Refs. 16, 
17] develop an extended Huygens-Fresnel theory by adding a 
random phase term to the Huygens-Fresnel integral. In fact, 
it is equivalent to applying the Rytov approximation to the 
integral equation. After applying the paraxial approxima- 


tion, the extended integral equation becomes 





: (ik|r-r'|) = 
E(r) = 5 | Se oa D) 
S A 
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where y is a complex phase (equivalent to S1: Yura's 
notation), the primes denote aperture coordinates and the 
non-primes denote image plane coordinates. The Huygens- 
Fresnel integral is developed in standard texts like 
[Ref. 18]. We should note that (2.42) reduces to the 
Green's function solution of Tatarski and Clifford in the 
geometrical optics limit (i.e., v = k fn, (z)dz). This is 
not surprising Since the kernal of the Huygens-Fresnel 
integral is a Green's function. 


Yura now finds the average intensity <I(r)>, 


51044... ...5 
amc ff ә (EE EIS 


|r-r"| 


x QU d^r'a?r" » . (0 
The only functions in (2.43) that depend on time are 
the complex phases v' and bn, Hence, evaluating the 
ensemble average, if we assume ergodicity, becomes an effort 
in evaluating the average of the v's. We now change 
variables to those of Yura [Ref. 17] and Walters [Ref. 14], 


+ — —= HA — — 
p = (r'+r")/2 , 6 zri MEM RK R: 


With this change of variables, the Fresnel approximation 


becomes 
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ә = a ste — ce — 
ит — ri je , 02944) 


= 
D 
(D 
X 
(D 
N 
| 
N 
z 
II 


0. This gives a very physical result for 


the average intensity 


=1 ko 





2 
aao) = a faro rare F (2.45) 
where 
– - (H 5:5) 
му (57) = fegt + poetet - je ? gei (2.46) 
and 
M_(57,2) = <elv'(F'F)+Y"(F",T)], , (2.47) 


Since 2 >> r' or r", the exponential in (2.46) can be 
approximated as unity. This makes Mr the autocovariance of 
the aperture (or if normalized, the aperture mutual 
coherence function (MCF)). ME is the atmospheric mutual 
coherence function. This allows us to treat the propagation 
problem in two distinct parts; the effects of the atmosphere 
and the effects of the receiver optics. 

The trick now is to find the atmospheric mutual 
coherence function. Yura assumes that the amplitude 
fluctuation, x, and phase fluctuation, s, that compose v are 


Gaussian variables, which implies that exp(¥'+¥") is log- 
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normally distributed. Using a property of log-normally 
distributed random variables and some results of Fried's 
work in [Ref. 19], Yura represents the atmospheric MCF in 


terms of the wave structure function D(p), giving 
1 "x 
o VEREIS. =D MS (2.48) 


where D(p) = D, (0) t Од (о). Now we can apply the results of 
Tatarski's work, since we have forms for Р, (К,0) and ЕР. (K,0) 
which can be related back to their structure functions. 


From this, D(p) becomes 


L 
D(p) = dai x C SL sp UND (2.49) 
0 n 0 
0 
and 
L 
D(p) = 2555 for dz s POŻAL $ (2. 59) 
0 


where 0 is the inner scale and VAL is the Fresnel zone 


Size. Based upon (2.49) and (2.50), Yura represents (2.48) 


as 


(903977), = MCF(p) = exp (55) 2/3 , (2.51) 


0 


Se 
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where PO is the transverse distance for the MCF to decline 
by el, yura has developed forms for py including inner and 
outer scale effects [Ref. 17]. We should remember that PO 
is a distance in the aperture of the optical system. It is 
nominally on the order of a few centimeters during the day. 
4. The Modulation Transfer Function Approach 

Fried in (Refs. 19, 20] uses the modulation transfer 
function (MTF) rather than the MCF. This application of the 
MTF and MCF are closely related since both represent the 
effects of the atmosphere in different planes (i.e., 
aperture and image planes). The MTF is the modulus of the 
complex optical transfer function from the linear systems 
approach to optics. This points to the fundamental 
difference between Fried and Yura's approach.  Fried's work 
is done in the coordinates of the image plane, i.e., spatial 
frequencies, while Yura works in the aperture. Both are 
equivalent by the Wiener-Khintchine theorem since a lens 
produces the Fourier transform of the incident electric 


field in the image plane. The transformation between the 


two planes is 
p> ARK , 


where Res the focal length of the optics and K is thej spatial 
frequency of the scintillations. This allows us to move 


between the aperture and image planes. Thus, the MTF(K) is 
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RR 3 


MTF(K) = exp[-3.44 |, (29527) 


where ro is the 1/e distance of the MTF and is given by 


Lo = (ggg) k* sec(%) oss à; (2.53) 
0 


where ¢ is the zenith angle [Ref. 21]. 

Explicitly in the MTF and implied in the MCF is a low 
pass filtering of the scintillation spatial frequencies by 
the atmosphere and optics. This linear systems view is 
extremely powerful. It allows us to decompose the turbulent 
propagation problem into a product of atmospheric and 
optical system filtering functions. We will use this idea 
and the Wiener-Khintchine theorem later to our advantage. 

5. Measures of Spatial Coherence 

As mentioned in the introduction, PO and ro 
represent two measures of the atmosphere's ability to 
maintain the spatial coherence of a propagating wave. 
Nominal values for rg are approximately 3-10 cm for ground 
to space measurements. From (2.53), we see that turbulence 
close to the aperture of an optical system contributes most 
to rg Since Cə a is largest in the boundary layer [Ref. 2]. 

Fried in [Ref. 21] derives the theoretical equation 


for another measure of spatial coherence in the atmosphere. 
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This is the isoplanatic angle, d Isoplanatism is the 
dependence of the optical transfer function of a system on 
the position of the source. Since, Yura and Fried's 
approaches break the propagation problem into MCF's and 
MTF's for the atmosphere and optics, the atmospheric 
turbulence does not affect the OTF of the optical system. 
It affects the OTF of the atmospheric part, hence the term 
anisoplanatism. This means the atmospheric OTF is dependent 
on the particular propagation path through the turbulence. 
The isoplanatic angle is an angular measure of 
spatial coherence. If we consider two paths through 
turbulence in the atmosphere, the isoplanatic angle relates 
the mutual coherence . point betveen the tvo paths. Fried 
defines the function S(T,6) given by 


2 5273 


225 — 2.905 k ave ?(irti- (201777 "EUN 


path 


- liogn-Q11? « 2eeni- QDdve * (v)219/6 - 3i cen- qn1)* 


- 2r0[1-(2)]ve + {0v} 6) , (2.54) 


where L is the path length, v is a position along the path 
and c is the cosine of the angle between r and 9. In the 


limit as r/9 + æ (or two rays crossing at infinity), 


Dal 


ion S(Y,8) = 2.905 k? ә avè Paroa: ‘ (2455) 


r/0 + co 
path 
Fried defines 80 as 
ə. - 2.905 k* f ave,” v3 , (2056) 


path 


Interestingly, Eo is given by the other limit of the 


function S(r,9), 


lim S(F,6) = 2.905 k* J ave,7trt = wy o 
g path 

where 
tg > = (gg) k* fac 2 T (2.58) 


path 


and the (1 - y/L) ?/ 3 term is a spherical wave weighting 
function. Clearly, we can change variables from v to z by 


correcting for zenith dependencies, i.e., 


V PRZ SQ EE DEM 


where $ is the zenith angle. Now one can immediately see 
the relationship between rg and 9g. The isoplanatic angle 
measurement looking up is the same as the rg or MIE 


measurement looking down from space. 
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Walters has shown 80 relates back to PO and the 


Huygens-Fresnel picture. The relationship is given by 
89 = pg/L : (2-59) 


Using the Huygens-Fresnel picture, one can gain a physical 
insight into ro and e measurements. Consider a Source on 
the earth propagating spherical waves into space. 8g is the 
angle such that the MCF of the field at two different points 
on the wavefront is et icsmoriginal value: This picture 
can be reversed where two points in space emit spherical 
waves that sum at a point on the earth. This gives the 
corresponding rg measurement. One of the advantages to the 
Huygens-Fresnel approach is this reciprocity. 

Based upon the isoplanatic angle weighting function, 
turbulence between 8-15 km is weighted heavily. Hence, 
turbulence near the tropopause tends to make the isoplanatic 
angle Smamie (1.€., 3-5 urad). Large isoplanatic angles are 
in the neighborhood of 13-20 urad. The largest isoplanatic 
angles (20 urad) are very small, consequently adaptive 
optics concepts that compensate for the atmosphere have a 
monumental task since the correction is only good for an 
angle 95- 

6. Intensity Fluctuations 


Tatarski in [Ref. 6] obtains the variance of the 


intensity fluctuations from 


59 


5 -— , (2.20) 


vhere 7: is the covariance function at p = 0. He 
evaluates the integral over K in В, (0) for point sources and 
point receivers giving results for plane and spherical 
waves. In reality, we do have point sources, stars. 
However, real point receivers do not exist. Fried [Ref. 
22], obtains o” for finite apertures. He finds that the 


aperture averaged variance is 





2 2T D 
o > pł end MTF(p)B (o) , (2180115) 
X 
or 
4 0 


vhere D is the diameter of the optics and o is a constant 
depending on the type of aperture. — can be represented 
by the tvo-dimensional Hankel transform of the spectrum 
because isotropy imposes cyclindrical symmetry. 


Substituting this into (2.61) gives 


D 
>)" Ј ədə MTF(0) [2m TI iepel (K,0)KdK]. (2.62) 
TD 0 0 X 


2 T 


g = 





( 


si 


The p integration is trivial since we have the Hankel 


transform of an MTF or autocorrelation. This transform 
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yields the intensity spectrum for the finite aperture 


involved, I(K). Hence, the normalized variance is 








9.” ҝ 
5 — P f KdKF (K,0)I(K) ; 1225 
S X 
0 
where 
s2 = en; 2 
> 4 0 P 


Again we see the filtering of spatial frequencies by 
the atmosphere in P (K,0) and by the optics in I(K). 
Aperture averaging is the filtering done by the optical 
System. Equation (2.63) has limits for point apertures and 
infinite apertures. In [Ref. 6], Tatarski shows that for 
point apertures, (2.63) behaves as 27/6, while for infinite 
apertures it behaves as 29, This implies that all other 
apertures fall between these two z dependencies. 

Over the years, research has modified and shown the 
limits of Fried's initial aperture averaging work [Ref. 22]. 
Yura [Ref. 24] suggests a form for the covariance function 
for arbitrary z. This modified Fried's results slightly. 
However, a more serious problem was found by Dunphy and Kerr 
[Ref. 25] and by Homstad, et. al [Ref. 26], if turbulence 1s 


very strong. Experimental work has failed to show the 


theoretical aperture averaging effects in strong turbulence. 


41 


The trend is towards less aperture averaging than expected 
in strong turbulence. Azar, et. al. [Ref. 27] and Clifford, 
et. al. [Ref. 28] explain that the covariance function for 
Strong turbulence is governed by two transverse scale sizes. 
The short size determines the initial fast drop in the 
covariance function, while a second scale characterizes the 
long tail. The long tail of the covariance function reduces 
the effects of aperture averaging. The aperture averaging 
controversy will arise again when we look at isoplanometer 
data. 
7. Problems in thealheory 

The theory for wave propagation presented so far has 
been for weak turbulence or single scattering. We have 
neglected the effects of multiple scattering. When the path 
lengths are long or the turbulence is strong, the linear 
theory breaks down. Experiments show that for intensity 
variances of 0° ~ .3, the normalized variance saturates to l 
[Ref. 28]. Multiple scattering destroys the spatial 
coherence of the wave. The linear theory assumes a coherent 
wave incident on the perturber whose strength is dependent 
on ода (or the temperature gradient). If the spatial 
coherence is degraded such that it is less than the Fresnel 
zone size, VAL , the wave is partially coherent across the 
perturber. This implies that interference can no longer 


completely modulate the wave. Hence, the variance of the 
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fluctuation intensity tends to saturate or approach a 
constant. Other theories have evolved to handle propagation 
in strong turbulence or multiple scattering. Many of the 
theories find their basis in techniques developed for 
quantum mechanics. Stroebehn in [Ref. 29] reviews iterative 
Green's function solutions that give an exact answer to the 
multiple scattering problem. But the solution, in practice, 
cannot be computed. He also summarizes the concept of 
Feynman diagrams originally developed for quantum electro- 
dynamics, and other more sophisticated theories developed 


for strong turbulence. 
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III. REMOTE MEASUREMENTS OF THE ISOPLANATIC ANGLE 


A. TECHNIQUE 

Loos and Hogge [Ref. 23] extend Fried's work of [Ref. 
22] by completing the Hankel transform of the MTF for a "top 
hat" aperture function. Loos and Hogge [Ref. 23] and 
Walters [Ref. 30] realized that by measuring the normalized 
variance of the intensity fluctuations of stars one could 
obtain the isoplanatic angle to within a constant. Walter's 
technique is as follows. The normalized variance of 


intensity fluctuations is 





2 L A 
2 
G 
> = 4(4r°)(.033k°) f e, ma J kako(k)sin? k) , 
S 
0 0 


(Snr) 


where K is the spatial frequency of the amplitude 
scintillations, L is the path length over which or is non- 
zero, I(K) is the intensity spectrum of the receiver optics 


and (K) is the xı 


portion of the Kolmogorov spectrum 
for turbulence in the inertial subrange. The similarity to 


the isoplanatic angle given by 


L 
=- =. 2; 905 x [sec иә f 226,2 (0 227? MECZ) 
0 
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is obvious. The trick is to make the weighting function 


v(z), 
` 2 K”z 
W(z) = J xaxo(x) sin (RI (K) ; 57. 
0 
5 š e 
behave as z . Mathematically this means that 
2213 ei АМ(2) , (3.4) 


where А 15 a proportionality constant. Hence, the 


isoplanatic angle becomes 
-—-. 5 


vhere A" is a nev proportionality constant taking account 


mür and Pi's, etc. of (3.1). 
Previously we noted that the weighting function integral 


5/6 


behaves as z” for infinite apertures and z for point 


apertures [Ref. 6]. Luckily, the isoplanatic angle 
weighting function, 25/3. lies within these two limits. 
Hence, theoretically it should be possible to mimic this 
behavior. 


The power of this isoplanatic angle measurement tech- 


2 : 3 
nique is that no explicit knowledge of the Ca profile rs 
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required. Ca in the free atmosphere is difficult to 
measure remotely with high vertical resolution (as we shall 
see later). Hence, integrated parameters like ro and 80 
depending on cha are usually measured independent of a ae 
profile. Another important advantage that naturally evolves 
from this technique is the filtering done by the aperture 
averaging. The averaging aids in suppressing the high 
Spatial frequency scintillations that are most affected by 


passage through the atmosphere and are most liable to cause 


Saturation. 


B. THE FIRST ISOPLANOMETERS 

Walters developed the first isoplanometers to operate 
during the day and night [Ref. 30]. The optical receivers 
an Celestron 14 Schmidt-Cassegrain telescope obscured so 
that a 11 cm diameter clear aperture remains. The theory at 
the time of instrument development was based on a "top hat" 
aperture function. Walters found that the 11 cm aperture 


gave a good approximation to the z2/3 


weighting function. 
Since these instruments are designed to operate at both 
day and night, they have a very small field of view to 
suppress the Rayleigh scattered background encountered in 
daytime viewing. Hence, high quality telescope drives and 
mounts are required for accurate tracking. In short, the 
first three generations of isoplanometers require full-time 
support. In [Ref. 30], Walters shows that these instruments 


operate properly within the weak turbulence limits. 


46 


C. APERTURE APODIZATION FOR WEIGHTING FUNCTION IMPROVEMENT 

1. Weighting Functions for the Top Hat Aperture 

One can compute numerically the integral over K in 

(3.3) giving the theoretical weighting function for a 
certain I(K). In reality, the limits to the K integral are 
not zero and infinity. The spatial frequencies of the 
Scintillations are truncated by the inner and outer scales 
of turbulence. Hence, the upper and lower limits of 


integration become 


ur. 
Kmax t 
and 
7: 
Knin ^ Lo ; 


where LO is the inner scale and Lg is the outer scale. For 
the following computer simulations, ły - 5 mm and Lg = 10 m. 
The integration over K is done by a brute force application 
of Simpson's rule over 512 points for a given value of z. 
The z stepsize is log (Az) = .05. 

Figure 3.1 shows the error from - for a 10 cm top 


hat aperture as a function of z for a 500 nm electromagnetic 


vave. 1(K) for the 10 cm top hat is the Airy pattern, 


20, 823 | * 
I(K) = —KD. . , (250) 
(> 
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were sD is the drámeter of the receiver optics. “The curve 
in Figure 3.1 is normalized so that the integrated error 
between 2 and 20 km is zero. This ensures that the error is 


DA 


not consistently above or below z giving erroneous 
results when this weighting function is integrated into a 
measurement system. The normalization interval 2 to 20 km 
results from the fact that turbulence, hence can is large 
enough in that region to contribute to the isoplanatic angle 


pA: dependence. Clearly, Coe is largest 


integral given the z 
in the boundary layer [Ref. 2], however the weighting is 
very small there causing little change in the isoplanatic 
angle. Actual profiles of Gu like that provided by Dr. E. 
Good from the Air Force Geophysics Laboratory, Hanscom AFB, 
MA in Figure 3.2, bear out this assumption. The particular 
profile shown will be used extensively in Chapters IV and 
VIII of this dissertation. A more complete description of 
the profile will be provided at that time. 

En [Ref. 24], Loos and Hogge plot the two 
theoretical limits for point and infinite apertures as a 
fiumetion of.z. They also plot the aperture averaged 
weighting function for a 35 cm telescope. They note that 


> In Figure 3.4 


2 
we see two different slopes present. If W(z) = z , the 


5/6 
slope would be +1/3, while if W(z) = 2 / , the slope would 


: 2 : 
the large aperture is on the z Side of z 


22775510O At low z, the slope in Figure 3.1 is 0s oo IV DP 
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at high z the slope is -0.63. As the diameter gets larger 
the slope of the low z portion moves closer to 1/3. 
Clearly, these figures only show the general influence of 
the limits. Figure 3.3 shows the movement of the peak 
towards higher z with increasing aperture size. Thus, we 
see more z” behavior as expected. 

The question arises as to why, with increasing 
aperture size, one sees more z^ behavior. The answer is 
straightforward. From [Ref. 23], we see that as D, the 
diameter of the optics goes to infinity, I(K) = 6(K-0). As 
finite D increases, the width of the Airy disk narrows due 
to the Hankel transform. Eventually, the Airy disk becomes 
a good approximation to the 6-function, hence, the z^ 
behavior. 

Figure 3.4 shows the wavelength dependence of the 
weighting function. To first order, the normalized variance 


given by (3.1) is wavelength independent. This can be shown 


by expanding sin^(K^z/2k) in a Maclaurin series giving 


£(z) = sin’ (=) X om) 55. (3.7) 


lle k? cancel in (3.1) removing the wavelength dependence. 
2 : : | 
We also see from (3.7) how the z weighting function 


mentioned above comes about. The weighting function, (3.3), 
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is itself not wavelength independent, hence we see the 
behavior persent in Figure 3.4. As the wavelength increases 
the weighting function moves laterally towards lower Z. 

The calibration constant, A", from (3.5), for the 
10 cm top hat aperture is nominally 1.0 x 1010, We would 
expect that this constant should get larger as more aperture 


averaging is included since the normalized variance will 


decrease. 


2. Improved Weighting Functions 
5/3 


The error from z in Figure 3.1 at 100 m is about 
a factor of three. In the weighting function given by 
(3.3), the only variable that depends on the measurement 
system is I(K). The remainder is due to the amplitude 
Fluctuation spectra F,(K,0). For some time, different 
groups have been modifying telescope apertures or image 
planes to measure GNE profiles [Refs. 31, 32, 33]. The 
modification of aperture functions by masks, etc., is called 


5⁄3 weighting function 


aperture apodization. Since the z 
lies between the theoretical limits, we hope that through 
aperture apodization we might improve the isoplanometer 
weighting function. 

To obtain intensity spectra from different aperture 
functions, one can Hankel transform the MTF or autocorrel- 


ation of the aperture or Hankel transform the aperture 


function and square the resultant amplitude spectra. The 
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results are equivalent by the Wiener-Khintchine theorem. 
Bateman [Ref. 34] compiled many different Hankel transform 
pairs. Using this list and different combinations of Bessel 
functions, I began constructing different weighting 
functions. Appendix A is a synopsis of the some apertures, 
spectra and resultant weighting functions. We found that 
combinations of Bessel functions are most desirable for 
spectra, since they represent apertures that can be 
implemented easily on a telescope. 


Before moving on to the aperture function that 
5/3 


5 


approximates the z quite vell, it is interesting to look 


asthe error from z produced by a simple annulus 
aperture. Figure 3.5 is an error plot for an annulus vith a 
10 cm outer diameter and a 4 cm obscuration diameter. The 
tvo peaks correspond to the influence of the inner (lover Z 
peak) and outer (higher z peak) diameters of the tvo Bessel 
functions involved. This leads one to believe that by using 


a dependence could 


combinations of Bessel functions the z 
be mimicked. Figure 3.6 is the error for a Celestron C8 
telescope aperture. It is interesting to note that the 
inner diameter affects the weighting function at high z, 
while the outer diameter affects performance at low z. We 
saw this before with the two theoretical limits of the 
weighting functions. 


After some effort, an aperture function was 


identified for a 8" Schmidt-Cassegrain telescope that 
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produces a weighting function that approximates the zd 


dependence. The aperture is the double annulus shown in 


5/3 


Figure 3.7 producing the error from z shown in Figure 3.8 
for a 500 nm electromagnetic wave. Over the altitude range 
critical for isoplanatic angle measurements (2-20 km), this 
weighting function does very well. 


The weighting function for the double annulus 


aperture function is given by 

















2 e KD 
1 f m. 
W(z) = e ss Fan sins (512 A 25) 227 
IAS ZET =. 0 
2 
e,KD e Ë KD 2 
w Gs 2:Ф6| |. z] Wz 3.8 
1 e j KD KD 3 c KD , (3.8) 
2 : 2 


where €), £5, €4 are constants that ratio the inner 
diameters to D, the outer diameter of the optics. (K) is 


the ka. 


dependence of the Kolmogorov spectrum. Figure 
3.9 shovs a comparison of the amplitude spectra squared 
(intensity spectra) for the 20 cm double annulus and a 20 cm 
top hat aperture. It is interesting to note that the 
increased higher spatial frequencies present in the double 


annulus spectrum are responsible for the improvement in the 


weighting function. 
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3. Double Annulus Calibration 
The calibration constant, A', for the 20 cm double 

annulus aperture function is approximately a factor of four 
larger than the 10 cm top hat. This is due to the increased 
aperture averaging. The double annulus weighting function 
shows the same wavelength dependence as the top hat, Figure 
3.10. An increase in the wavelength of the EM radiation 
shifts the weighting function to lower z. One interesting 
point is that the normalized variance has no wavelength 
dependence as we have shown. However, the isoplanatic angle 
certainly does (3.2). Clearly, all the wavelength 
dependence for the isoplanatic angle is in A' and it should 
behave, to first order as k2. Figure 3.11 shows the A' 
wavelength dependence. At wavelengths greater than 800 nm, 


R but diverges belov 800 nm. Clearly, this 


A" behaves as k 
is important depending on the temperature of the stars used 
for isoplanatic angle measurement. 

Based upon the factor of four increase in the 
calibration constant, one might expect that an instrument 
vith this veighting function vould be able to vithstand a 
factor of four more turbulent intensity before saturating. 
However, [Refs. 25, 26, 27] show that the aperture averaging 
effects may be degraded in strong turbulence making the 10 
cm and 20 cm devices saturate at the same relative 
normalized variance. This issue will be resolved later by 


experiment. 
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The 20 cm double annulus has approximately twice the 
areamol she 10 cm top hat. This implies that the telescope 
will collect twice the photons increasing the average 
intensity by a factor of two. The variance of the shot 
noise current will also increase by a factor of two. 
Mathematically, the isoplanatic angle for the 10 cm system 


1s 


where e," is the variance of the signal and e, is the 
variance of the noise. For the 20 cm double annulus system, 


this same equation becomes 


53/3 _ : 
0 = “Ar 0 


S 4I 


The proportionality constant, A', for the 20 cm system is 
approximately a factor of four larger than for the 10 cm 
system. In practice, we measure the shot noise experi- 
mentally and remove it from the variance in software. 

Hence, its effect is negligible at night. The 10 cm 
isoplanometer uses a pellicle to split the light beam in the 
detector subassembly. This reduces the average intensity by 
a factor of two. The 20 cm system, without the pellicle, is 


equivalent to the 10 cm system. 
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4. Double Annulus Improvements 

The single greatest advantage the double annulus 
aperture function isoplanometer would have over the previous 
Systems is that the 35 cm telescope is no longer required. 
This fact alone means saving approximately $8000.00 per 
System, not to mention that an isoplanometer based on an 8" 
telescope can be set up and operated by a single individual. 
The weighting function itself is a significant improvement 
over that of the top hat. We will see in the next chapter 
just what this improvement means in terms of isoplanatic 


angle measurement. 
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IV. APODIZED ISOPLANOMETER SENSITIVITY ANALYSIS 
A A тә AO LS 


Two questions might be asked based on the double annulus 
PEEELMETNGP Funeeion ef Figure 3.8. First, at 100 m there 
remains a factor of 1.75 error. e at low altitudes is 


348 


large, So deviations from z may cause significant errors 
in the isoplanatic angle measurement. Second, in the 
numerical calculations of the weighting functions, a fixed 
inner scale was used for all z. Clearly, the inner scale is 
a function of altitude which may drastically change the 


weighting function. These two issues are addressed in this 


chapter. 


A.  WEIGHTING FUNCTION PERFORMANCE IN STRONG TURBULENCE 
Consider an isoplanometer viewing a star directly 
through a convective plume. Clearly, the values of CN will 


be very large; on the order of 10712 m 2/7 


near the ground. 
From Figure 3.8, one can immediately see that the weighting 
function for the double annulus has its greatest error in 


3/3 weighting is low near the 


this region. Even though the z 
earth, the = values in strong turbulence may be large 
enough to cause significant error in isoplanatic angle 
measurements. Based on this concern, I performed a 


sensitivity analysis on both the 10 cm top hat and double 


annulus aperture weighting functions. 
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1. The n Profile 


The most realistic way to test the weighting 


functions is to use real ә data. Figure 3.2 shows Good's 


Ce data between 500 m and 40 km. Good obtained this 
profile by differential temperature measurement between two 
balloon-borne microthermal probes. From these differences 
one can construct the temperature structure function DT and 


2 2 


obtain or as a function of z. ES and CT are related by 


the expression, [Ref. 6], 


AT 


ci = с әрли х 1079)2, 5— 


vhere P is the pressure in millibars and T is the 
temperature in Kelvin. Direct measurement of e in this 
manner yields very high vertical resolution, on the order of 
tens of meters.  AFGL averaged raw data to give the 500 m 
resolution in Figure 3.2. 

For purposes of this analysis, the region of 
greatest interest lies below 500 m. Walters [Ref. 2], plots 
c versus altitude for 1 to 1000 m. Walters measured this 
particular data set at midday in a desert location and it 
represents very strong turbulence conditions. To account 
for convective plumes in low altitude turbulence, I doubled 


2 


the Walters ską values. This gives one Ch profile from 1 m 


to 40 km. 
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2. The Analysis 


For convenience, I segmented the profile into four 
regions, l to 10 m, 10 to 100 m, 100 to 500 m, and 500 m to 


40 km. Over each of these regions, I numerically evaluated 


IE 2 (2) z?/ 3az P 
n 


2 
fe, (Z) Wny (Z) dz ; 


VESTE e 


Based on Figures 3.1 and 3.8, we would expect the weighting 


three integrals, 


Theoretical 


(4.2) 


Top Hat Aperture 


Double Annulus Aperture 


function integrals to be smaller than theoretical at low z. 
Table 1 shows the results of this analysis. Table 2 
contains the relative errors for each segment and the total 
path. The integrated values in Table 1 show that even the 
strongest turbulence at low z accounts for only a small 
percentage (about 2%) of the integrated path effect. 
Clearly, the errors over the total path for both the top hat 
aperture and the double annulus are within the limits of the 
theory. As one might expect, this analysis is extremely 
Sensitive to the choice of A, a proportionality constant, 
where 


ә 251 (4.3) 
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TABLE 1 


CONTRIBUTIONS TO THE ISOPLANATIC ANGLE INTEGRAL FOR 


Theoretical Top Hat Double Annulus 
1-10 m o вкл әса 5.80 x 10 7” 
10-100 m 8.19 x Low? E TE ` una 
100-500 m I IGE s. 50 ә 8.67 x 107” 
500 m-40 km 3.89 x 107° ocio 3.91 x 1077 
TABLE 2 
RELATIVE ERROR FROM THEORETICAL FOR STRONG 


= LC jg 
10-100 m 
100-500 m 
500 m-40 km 


1 m-40 km 


TURBULENCE ANALYSIS 


Top Hat 


Fo 
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SƏTRONG TURBULENCE ANALYSIS 


Double Annulus 
. 78 
.54 
226 

6 x 10. 


SW TO 


In building an isoplanometer, the better one can define A or 
A', the better the instrument should operate. 

This analysis shows that an isoplanometer should 
perform well, based on its weighting function in strong 
turbulence. However, we must note that in very Strong 
turbulence, saturation has already reduced the reliability 


of the isoplanatic angle measurements. 


B. INNER SCALE EFFECTS 


The inner scale, n, is given by 


n = (>) (4.4) 


where v is the kinematic viscosity and e is the energy 
dissipation rate. Clearly, the kinematic viscosity and the 
dissipation rate are both functions of altitude. In com- 
pleting the weighting function integral (3.3) numerically, I 
used the same inner scale for all values of z. This could 
mean that the weighting function does not accurately model 


the 2713 


behavior required for isoplanatic angle 
measurement. 


l.  Constructing the Inner Scale as a Function of 
Altitude 


From (4.4), we see that the inner scale depends 


strongly on the kinematic viscosity v, given by 


del 


where 6 is the viscosity coefficient and p is mass density. 
The U.S. Standard Atmosphere [Ref. 35] has vertical profiles 
of the kinematic viscosity. Hence, this portion of the 
inner scale is well understood. 

The dissipation rate is another story entirely. 
Turbulence in the free atmosphere is intermittent, 
temporally and spatially, unlike boundary layer turbulence. 
This is compounded by the fact that measurements in the free 
atmosphere are difficult to make. Using [Refs. 36-41], I 
constructed what might be called a "nominal" dissipation 
rate profile. This profile, by no means, represents the 
output of any theoretical model on my part. It 1S a 
synopsis of the data and model results given in [Refs. 
36-41]. The "nominal" dissipation rate profile is shown in 
Figure 4.1. Some controversy exists as to the behavior of 
the dissipation rate at higher altitudes [Ref. 38]. It is 
generally agreed that e can reach values comparable to that 
observed in the surface layer. My profile shows light- 


3 (or 


moderate turbulence near the surface, e = 40 cm^/sec 
log e = 1.6) at 100 m and decreases to .38 en Sec (or 

log € = -0.42) at 55 km. As a point of reference, we might 
note that the dissipation rates in convective storms are on 


3 4 


the order of 10-10 [Ref. 40]. The units on 
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dissipation rate indicate a "specific" dissipation rate 
where system mass is removed. 

Based on the Standard Atmosphere kinematic viscosity 
profile and my "nominal" dissipation rate profile, one can 
now construct the behavior of the inner scale with altitude. 
Figure 4.2 shows the inner scale profile. Clearly, the 
inner scale changes dramatically with altitude. At 55 km, 

n is on the order of a meter. As expected, the inner scale 

follows the kinematic viscosity since kinematic viscosity 

increases three orders of magnitude between 100 m to 50 km. 
2.  WeightfIng Functlon DependencemonmthcasmH NE SN 

A crude approximation to the inner scale effect can 
be obtained by increasing the inner scale.for all z while 
numerically evaluating the weighting function integral. 
Figure 4.3a shows how the double annulus aperture weighting 
function changes with an inner scale of 10 cm. The 
deviation in W(z) caused by the change in inner scale from 5 
mm to 10 cm is not catastrophic, however, a more complete 
analysis needs to be done. 

Since the upper limit of the K integration depends 
on the inner scale, we see that increasing inner scale size 
corresponds to lowering the upper limit of integration. We 
expect that W(z) will not change significantly at lower z, 
due to the modest increase in the inner scale. However, as 


the inner scale increases, the upper limit of integration 
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will decrease. At some critical K the upper limit will 


max?” 
begin eliminating integrand producing a dramatic change in 
the weighting function. 


5/3 for the double 


Figure 4.3b shows the error from z 
annulus aperture weighting function with the z dependence in 
the upper limit of integration. As expected, we see little 
change at lower altitudes. At 40 km, the inner scale is now 
large enough to begin removing part of the integrand. The 


5/3 


deviation from z begins when the inner scale reaches 15 


5/3 is a factor of three. 


cm, and by 25 cm, the error from z 
The inner scale behavior of the weighting function 
is not a problem in a measurement system, since values of 
SE at 40 km are. approximately NUS nü Disturbances 
to the electromagnetic vave propagating through the medium 


2: 
occur at lower z where C. 1s larger. 


C. SENSITIVITY ANALYSIS CONCLUSIONS 

Based on the results of sections A and B, the 
isoplanatic angle weighting function appears very robust. 
The only real theoretical problem, in terms of measurement, 
is saturation of the normalized variance. One might ask 
about outer scale effects.  Nominally, the outer scale, 
except very near the surface, remains on the order of 10 m 
[Ref. 42]. This can be seen from the Reynolds number. The 


Reynolds number is, (as shown in the Background section), 
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mere v is the velocity, L is a characteristic size of the 
flow and v is the kinematic viscosity. The onset of 
turbulence occurs at Reynolds numbers of 2000 or greater. 
in the free atmosphere remains about the same and v 
increases with altitude. Thus, velocity is the parameter 
responsible for large Re. Since L does not vary 
considerably, Lo” the next scale size smaller, should not 
vary considerably either. Near the earth, L has an 


immediate fixed boundary condition that reduces its size. 


Tə 


V. INSTRUMENT QUALIFICATION 


Currently, many DoD agencies require isoplanatic angle 
measurements at night. Based on this need, we fabricated a 
fourth generation isoplanometer. What follows is a brief 
description of the hardware, detailed analysis of instrument 
qualification and a comparison with the 10 cm "top hat" 


isoplanometer. 


A. HARDWARE 
l. Telescope 

The double annulus aperture function of Figure 3.8, 
was developed for use on an 8" Schmidt-Cassegrain telescope. 
The Celestron C8 is a reasonably priced and readily 
available Schmidt-Cassegrain telescope with a high quality 
worm-gear drive for star tracking. The choice of Bessel 
functions as a spectral basis set has the advantage of 
simple obscuration implementation. The aperture obscuration 
is two disks made of flat black rubber adhered to the 
Celestron corrector plate with rubber cement. This allows 
removal at a later date without damage to the optics or 
coatings. With the exception of the obscuration, all of the 
equipment associated with the telescope (tripod, mount, 


etc.) is furnished by Celestron. 
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ż. JaBocector 

The detector subassembly is a modification of the 
already existing 10 cm "top hat" isoplanometer detector 
developed by Walters [Ref. 30]. Walters engineered and 
implemented the modifications to the fourth generation 
detector subsystem, associated signal conditioning hardware 
and software. The detector subassembly screws into the rear 
of the telescope. Incident radiation is directed over one 
of two paths. To acquire stars and focus the telescope, the 
light enters a wide angle eyepiece. Otherwise, the light 
passes through a pupil (on the order of millimeters), and is 
incident on a photomultiplier tube that senses the intensity 
fluctuations (stellar scintillations). The image is 
slightly defocused to remove the effects of any inhomogen- 
eities on the surface of the photomultiplier. The photo- 
multiplier is chosen to have good red response, which causes 
a slight improvement in the weighting function at low alti- 
tudes, Figure 3.10 and suppresses any Rayleigh backgound. 

We should note that the pupil in this instrument can 
be much larger than the previous isoplanometer, since it is 
intended for night use only. The 10 cm isoplanometers use a 
chopper to make reference background measurements necessary 
for daytime operation. This feature is not required on the 
fourth generation system for the same reason as above. A 
Schematic diagram of the detector subassembly is shown in 


Figure 5-1. 
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3. Signal Conditioning and Isoplanatic Angle Extraction 


The first fourth generation isoplanatic angle 
measurement system (A system), uses a Hewlett-Packard 
HP 3421 Data Acquisition unit to sample the output of a 
signal conditioning unit. The signal conditioning unit 
performs two functions. It low pass filters the incoming 
signal below 500 Hz and calculates the mean and rms values 
of the input voltage data from the photomultiplier. The 
Data Acquisition Unit, with an internal voltmeter, detects 
the voltage on each channel once/second and relays the 
information to a Hewlett-Packard 217 Computer. 

Once the data is in the computer, the normalized 
variance of the signal is calculated from the mean and rms 
values of the voltage. Obtaining the isoplanatic angle is 
then a simple matter of applying the constant A' to the 
normalized variance. The average intensity (in volts) and 
isoplanatic angle (in urads) are plotted on a graphics 
monitor. A real-time plot of the average intensity can be 
very helpful during data acquisition in determining the 
occurrence of clouds and monitoring stellar tracking errors. 
The isoplanatic angle is very robust to drifts in the 
average intensity since we normalize the variance. 

On later systems (B systems), the HP 3421 is 
replaced by an Infotek A-D converter that resides Wl thie ae 


HP computer. The computer samples the A-D converter so that 
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one isoplanatic angle measurement is made per second. The 
computer does not store all the isoplanatic angle data. 
Instead, the computer calculates ten second averages and 
stores the averages every 15 minutes. 

The fourth generation isoplanometer uses software 
that has evolved and become more "friendly" with each field 
experiment. Currently, the zenith angle dependence of the 


isoplanometer angle is removed within the main program. 


B. ZENITH ANGLE QUALIFICATION TESTS 

It is an interesting problem to know if your instrument 
is really measuring the isoplanatic angle. One could fly 
balloons and obtain a direct measurement of the Chat profile 


5/3 


and integrate with z weighting, however, Walters [Ref. 
30], has developed a much less costly and simpler technique. 
The isoplanatic angle is 


/3 E 


-5/3 
d A" (3) , (5.1) 


80 = [sec ML 


including the zenith angle dependence. For increasing path 
length (higher zenith angles), the normalized variance will 
increase since the light is encountering more integrated 
turbulence. Hence, a plot of E versus log (sec 6$) 
should have an 8/3 slope if the instrument is truly 
measuring the isoplanatic angle to within a multiplicative 


constant. 
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Walters used this technique to qualify the 10 cm "top 
hat" isoplanometer (really 11 cm) [Ref. 30]. The following 
is an application of the above technique to the fourth 
Generation system», but first a caution. Clearly, if the 
atmosphere is horizontally homogeneous (i.e., composed of 
onionskin layers), the technique should be effective. 
However, if the atmosphere is irratic and the turbulence is 
not uniform over the night sky, this technique will not be 
as reliable. 

We (D. L. Walters and myself) obtained qualification 
measurements for the fourth generation system at two 
different locations. The first data set is from McDonald 
Observatory in western Texas. The night was ideal for 
instrument qualification using the zenith angle technique. 
The isoplanatic angles were on the order of 15y rad. In 
Figure 5.2 each point represents one 10 sec average of the 
normalized variance, while each group of points is 400 
seconds of data on a given star. Figure 5.2 shows very good 
agreement with theory on this test. The least-squares fig 
has a slope of 2.60, which is not a significant error from 
8/3. The least-squares fit does not include the two data 
sets at high zenith angles since the normalized variance is 
saturating. 

The second set of data is from Albuquerque, New Mexico, 


acquired when atmospheric conditions were less than optimal. 
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During the daylight hours, large thunderstorm cells formed, 
@issm@pating at night into a layer of cirrus. About 0200, 
the cirrus cleared and we obtained the data set in Figure 
5.3. The least-squares fit includes the three observations 
at lower zenith angles. The slope of the least-squares fit 
is 2.48, less than the 8/3 theoretical. Based on the 
conditions that night, the results are very good. 

More interesting in Figure 5.3 is the behavior of the 
normalized variance at large zenith angles (last two data 
sets). We note that the normalized variance is no longer 
changing for these data sets. Clearly, the path is long 
enough to induce saturation in the normalized variance. For 
the fourth generation isoplanometer this occurs at 
normalized variances of about .06 to .07. The 10 cm 
isoplanometer saturates at normalized variances of .2 to .3. 
The aperture averaging has lowered the normalized variance 
between the two instruments the predicted amount. Remember, 
the calibration constant A' for the 10 cm isoplanometer is 


approximately 1 x To vhile A" for the fourth generation 


isoplanometer is 4 x 0. Hence, for a given isoplanatic 
angle, the normalized variance for the nev system vill be a 
factor of four lover than the 10 cm system. Thus, the 
double annulus system is saturating at an equivalent 


normalized variance, We still have not resolved whether the 


increased aperture averaging aids in forestalling 
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saturation. The fourth generation system saturates at its 
nominal "saturation" Кәләми: variance, however, it may 
take four times the turbulence to make it saturate. The 
fourth generation system seems to take reliable data at 
slightly larger zenith angles than the 10 cm system. 
Resolving this question requires simultaneous observations 
with both instruments on the same stars for a zenith angle 
check. This test has not been possible due to instrument 


availability. 


C. APODIZED AND 10 cm ISOPLANOMETER COMPARISON 

Even though simultaneous zenith angle tests were 
impossible in Albuquerque due to instrument availability and 
atmospheric conditions, we did obtain simultaneous data on 
Vega on August 25, 1985. The data is contained in Table 3. 
The mean isoplanatic angle for the 10 cm isoplanometer data 
is 6.7 + .17; while the double annulus averaged 6.8 + .15. 
This shows very good agreement between the two instruments 
and removes the ambiguity of the multiplicative constant in 
the zenith angle test. Table 3 also shows an unstable 
night. The isoplanatic angles range from 4.5 urad to 9.0 
urad over a one hour period. 
DE ISOPLANATIC ANGLE MEASUREMENTS WITH SUPPORTING 

METEOROLOGICAL DATA 

During September, 1985, the fourth generation system was 


in use on Haleakalea crater, Maui, HI. One night, our data 
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Time 


(Local Standard) 


2210 
2212 
2213 
2214 
2216 
2219 
2220 
2222 
2224 
2220 
2228 
2229 
2230 
2231 


2233.00 


2233.41 


2234 


2255 


2207 


22516 


TABLE 3 


10 cm (urad) 


20 


10 cm AND DOUBLE ANNULUS COMPARISON 


Double Annulus (urad) 


TABLE 3 (CONTINUED) 


Time (Local Standard) l0 cm (urad) Double Annulus (urad) 
2240 53 6.9 
2241 6.3 7.4 
2242 ous 255 
2243 nü 6.6 
2246 nə. 11 
2247 6.9 6.9 
2249 T: о 
2250 Eer 5,2 
2252 5— 6.4 
2253 6.0 6.0 
2255 6.4 55 
2256 4.5 4.7 
2257 | 4.9 5 
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collection overlapped a meteorological balloon flight. The 
balloon provided altitude, pressure, temperature, dewpoint 
depression, relative humidity, wind speed, and direction. 
Figure 5.4 shows the real-time computer output on 17 Sep 
1985 (Star 17=Vega). The solid line is the time history of 
the average intensity, while the points represent 1 second 
measurements of the isoplanatic angle. The spikes in the 
average intensity occurred when the telescope dome lights 
were turned on. Figure 5.4 shows a small isoplanatic angle. 
However, no zenith angle correction has been applied to this 
data. Figure 5.5 contains this 15 minute set of data after 
application of the zenith angle dependence. The isoplanatic 
angles are in fact fairly large, approximately 13 urad. 
These large isoplanatic angles indicate low turbulence 
between 8-15 km. Hence, the wind profile should show low 
wind speeds and small shears. The complete set of meteoro- 
lecce data is contained in Appendix B, the wind shear data 
is as follows. The balloon data indicates three shear 
layers (2 km, 14 km, 27 km). The wind speed at the 2 and 14 
km levels is approximately 4 knots, while at 27 km is 47 
knots. These wind speeds are very low and indicate that 
very little shear driven turbulence is present. Hence, the 
large isoplanatic angles agree with the observed wind 


profile data. 
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Figur 5.4 Real-Time Output of the Fourth Generation 
Isoplanometer at Maui, HI 
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E. QUALIFICATION CONCLUSIONS 

Based on the zenith angle tests, the simultaneous 
measurements with the 10 cm device and the measurements with 
supporting meteorological data, we conclude that the 
instrument is indeed measuring the isoplanatic angle 
correctly. Currently, two fourth generation instruments are 
in operation at different locations worldwide. This data is 
being archived at the Naval Postgraduate School by D. L. 


Walters. 
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VI. DATA COLLECTED IN MAUI 


As mentioned previously, we conducted a trip to Maui, HI 
supplying an isoplanometer to the AMOS facility. During 
this trip, weather prohibited acquisition of large amounts 
of data, however, we did obtain that data mentioned in 
Chapter V, along with what I will present here. 

In this section, we will see one of the greatest assets 
of the fourth generation isoplanometer. That is its ability 
to make isoplanatic angle measurements effortlessly each 
second. The following plots show trends in isoplanatic 
angle data that have not been observed before. 

Figures 6.1-6.7 are real-time computer output of each l 
second isoplanatic angle measurement. As in Figure 5.4, the 
solid line is the average stellar intensity and each point 
is an isoplanatic angle measurement. We might note that the 
average intensity is much more steady than in Figure 5.4. 
This is due to the lack of high thin cirrus clouds. Also, 
in Figure 6.2, the steps are caused by changing the voltage 
to the photomultiplier. Note that these changes do not 
affect the isoplanatic angle measurements. In Figure 6.7, 
the deviations in the average intensity are due to the star 
beginning to move out of the telescope field-of-view. 

In these data sets, we see isoplanatic angles ranging 


from 16 urad to 5 urad. This, of course, indicates a very 
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Figure 6.1 Real-Time Isoplanometer Output from Maui, HI 
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Figure 6.2 Real-Time Isoplanometer Output from Maui, HI, 
Showing Voltage Changes to the Photomultiplier 
and a Rapid Decrease in the Isoplanatic Angle 
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Figure 6.3 Real-Time Isoplanometer Output from Maui, HI, 
Showing a Rapid Decrease in the Isoplanatic 
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Figure 6.4 Real-Time Isoplanometer Output from Maui, HI, 
with Oscillations of 90 Seconds 
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Figure 6.5 Real-Time Output of the Isoplanometer in Maui, 
HI, Showing Oscillations with 1 Minute Periods 
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Figure 6.6 Real-Time Output of the Isoplanometer in Maui, 
HI, Showing Oscillations 
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Figure 6.7 Real-Time Isoplanometer Output from Maui, HI, 
Showing the Star Moving from the Telescope 
Field of View 
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unstable night. But, what is most interesting are the low 
frequency oscillations present in the data. Figure 6.4 has 
oscillations with periods on the order of a one minute, 
while Figure 6.5 shows oscillations on the order of three 
minutes. Unfortunately, no balloon measurements of meteor- 
ological data were performed on 20 September. Without 
supporting meteorological data, it is impossible to 
ascertain exactly the mechanism responsible for producing 
data like this. 

Based on the isoplanometer weighting function, it is 
possible to make some general comments on the location of 
the turbulence that could produce this type of oscillation. 
Since the isoplanometer is weighted 22/3, low altitude 
turbulence (boundary layer to 2 km) is not responsible for 
this behavior. Most likely, the turbulence producing the 
oscillations is located in the tropopause. At 10 km, the 


z2/ 3 


n m n ore values are large enough to produce 
the lower isoplanatic angles we see. Researchers in 
Albuquerque, NM, indicate that they observe oscillatory 
data, with the second fourth generation isoplanometer, when 
the jet stream is moving into the area. As the jet stream 
approaches, wind speeds in the tropopause begin to increase. 
Periodically, the gradient Richardson number (ratio of 


potential temperature gradient to velocity gradient), will 


fall below the critical Richardson number (about R; = 1/4), 
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and the flow will give way to shear driven turbulence. The 
jet stream is not close enough to supply the energy 
necessary to maintain the turbulence. Hence, after a short 
period of time, the flow becomes stable and the isoplanatic 
angles increase. Once the jet stream has completely entered 
the area, the isoplanatic angles no longer oscillate, but 
remain low due to constant turbulence. Thus, one can 
surmize that the oscillatory isoplanatic angles are caused 
by temporally intermittent turbulence. 

More research needs to be done in the area of 
isoplanatic angle measurement with supporting meteorological 
data. Once trends in isoplanatic angle data have been 
established for certain atmospheric phenomena, balloon data 
may no longer be necessary for an accurate picture of 


turbulence in the tropopause. 
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VII. REFRACTIVE TURBULENCE PROFILING 


In the discussion on the isoplanatic angle, we saw that 
by measuring the normalized variance, we avoided the unknown 
da profile. Clearly, the integrated parameters like Bo and 
d do not provide direct information on the distribution of 
turbulence with altitude. Many laser system design analyses 
require SA profiles. However, they are difficult to obtain 
by remote measurements. In Chapters IV and V, I presented a 
profile measured by microthermal probes on a balloon. This 
represents a direct measurement of SE? with very high 


vertical resolution (on the order of meters), along a single 


2 


vertical path through the atmosphere. Actively, Ch 


profiles are measured by acoustic sounders (up to approxi- 
mately 1 km) and pulsed Doppler radars (2-30 km) [Refs. 41, 
43]. Active techniques provide good vertical resolution, 
but the cost of such instruments is quite high. 

In this chapter and the next, we will look at profiling 
one passively by direct inversion of the amplitude 
scintillation covariance equation with application of the 
Tikhonov regularization technique. But before delving into 


a rigorous solution, we should make note of the other 


passive techniques for measuring ue profiles. 
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HS PREVIOUS WORK 
The time lagged covariance function of amplitude 


scintillations for a point source and two equal circular 


aperture receivers is given by [Ref. 44], 





B (p!) = Ark? IER ] k&'ak'e ciao tkt po - v(z')«|] 
0 0 
0 
D e ии 
x 3 LEJE 5 Com 
2k K'D f Ы 
2 


where о" is the separation between the detectors, v(z') is a 
wind profile, t is the time-lag, 9?(K') is the 0.033k' 11/3 
portion of the Kolmogorov spectrum and D is the diameter of 
the circular receiving apertures. Often this equation is 
written for point receivers. In that case the Ji (x) /x ” 1, 
For measuring aha profiles, we let the time-lag, t, equal 


zero, removing the unknown velocity profile. Hence, Bin) 


is in the form of a Fredholm equation of the first kind, 
E NS Joe netos , 5 


where 6 3 constant and Q(z",o") is the kernel given by 
the integral over K". Fredholm equations of the first kind 
are notorious for being ill-posed (i.e., errors or noise in 


the data make the solution ill-behaved). In fact, Fredholm 
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equations of the first kind may have no solution at all. 
Peskoff [Ref. 45] was the first to invert this equation for 
point apertures. Shen [Ref. 46] realized that the problem 
was ill-posed and used a least-squares technique to minimize 
error. As time progressed, other methods for handling the 
ill-posed problem appeared. Heneghan and Ishimaru [Ref. 47] 
proposed a statistical inversion technique to compensate for 
the large errors induced by inversion. However, one should 
note that none of these techniques evolved into an actual 
measurement system. People began to realize that if the 
weighting function, W(z'), (3.3), could be made to peak, C 


n 
could be measured directly without inversion using (3.1). 

1. Remote Probing with Apodized Apertures 

To make W(z') peak requires filtering in the 

aperture or image plane [Refs. 31-33]. NOAA [Ref. 32] 
fabricated a зан profiling instrument (stellar scintillo- 
meter) based upon the filtering technique. By using three 
optical weighting functions, it measures values of Gun at 
seven levels between 0-20 km. Unfortunately, the weighting 
functions, W(z'), are somewhat broad. The structure we 
observe in Figure 3.2 would not be as pronounced if this 
profile had been measured by the optical scintillometer. 
However, the scintillometer does follow the general trends 
of NE well. 


In Appendix A, I show many weighting functions, 


W(z), for different apertures and spectra. This was an 
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effort to find highly peaked weighting functions for an 
optical scintillometer type instrument. Unfortunately, the 
quest proved fruitless, but other interesting weighting 
functions were found. 

2. Binary Star Techniques 

Several French researchers [Refs. 48-50] have 

developed a technique for profiling ee using a spacio- 
angular correlation function of two stars in a telescope 
aperture. The maxima in the correlation function correspond 


to levels of turbulent layers through the relationship, 


where 0 is the location of a correlation maxima, 90 is the 
angular separation of the binary stars and - is the 
altitude of the turbulence corresponding to the maxima. The 
group built instrumentation to measure che profiles in this 
manner. They report vertical resolution on the order of 2 
km, which is very good. However, being dependent on binary 
stars is a significant hinderance, esecially when the cha 


profile vertical resolution depends on their angular 


separation. 


Bs INVERSION WITH TIKHONOV REGULARIZATION 
Jarem [Refs. 3, 4] implemented a regularization 
technique for handling the ill-posed direct inversion 


problem on the covariance function for point detectors. 


PO? 


This has been the first "new" attempt at cla profiling 
recently. The Tikhonov regularization amounts to filtering 
the high spatial frequencies of the inverted spectrum. 
Since aperture averaging also filters high spatial 
frequencies, it seems reasonable to include this effect in 
the covariance function to aid in regularization. In the 
following sections, I will invert the aperture averaged 
covariance function by Hankel transform and apply the 
Tikhonov technique. This inversion is a general approach 
developed by Peskoff [Ref. 45]. The regularization portion 
follows Tikhonov [Ref. 5] and Jarem [Refs. 3, 4]. 
1. Inverting the Fredholm Equation 

The covariance function of amplitude scintillations 

for a single point source (star) and two identical but 


arbitrary apertures at zero time lag is given by [Ref. 44], 


со 





3 “ölə 
B (p!) = ant! | az'c_*(z') f "dk" 9(K")sin PZ) I (K') 
0 0 
X Jg(K'o') y (7:57 


vhere k is the vavenumber of the monochromatic electro- 


magnetic vave, K" is the spatial frequency of the amplitude 


27: 


scintillations, 9?(K') is the K^ dependence of the 


SITO 


Kolmogorov spectrum, p' is the detector separation and I(K') 
is the intensity spectrum of the receiver optics. We might 
note that the equation is valid only in the inertial 
subrange due to the Kolmogorov spectrum. Like Jarem, we 


will extend the limits of integration over = such that 
2 1 = 1 1 
C. EAS CL, 


where L is the path length over which ane is non zero. This 
equates to making Gren an even function. 
The two-dimensional Hankel transform of B_ (0°) is 


the spectrum z. i given by 


F (K') = 2“ E . (7.4) 
0 


EMbStMEUEiIng (7.3) into (7.4) gives 





23 2 | 2 | 2 K'^z 
EE) = 4T k fazre, (z') | aks (KTK) sine aK ] 
mares 0 
x İ ədə KO) (Kə ")90 (K"o") i (025) 
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Peskoff points out that 
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J p'dp'KJg (Ko')Jg (K'o') = 6(K-K") : (7.6) 
0 


This allow us to complete the K integration giving 


p - Av)k^I (K' ) 6 (K') Г дам E sin 


== СО) 





C. 


Using the double angle relation sin”a = (1/2) (l=cos2a), 


(7.7) becomes 


F_(K') = 2r3kÉ1(K')9(K') Ј ас̧а" a Paren a 


zz OO 





X COS [Az ] š (750 


Z 2 


Jarem notes that J aztec, (z') = 2c , where cha is the 


= me 


average value of 692 along the path. The second term 


represents the Fourier transform of c ^ (z'), which we will 


call C UQ substituting these results into (7.8) gives 


3 


E = 2 ' ' ız 
P UN ) = 2r k I(K')9*(K ) [2C 


2 (gt 
--— e (7.9) 


Jarem defines a set of dimensionless variables for 


convenience, 
2 ЈЕ 
x = K' (z) 3 uà 2 ; p = o 


where d is the maximum separation of the detectors. It is 
interesting to note that (L/k) is related to the Fresnel 


zone size since 


L 1/2 _ 


| REN: 
k 


( ' 
(217) 14? 


DF 


where f is the Fresnel zone size. With the new variables, 
$(K') becomes 
117: k -11/6 .-11/6 k, -11/6 


$(K') = .033K' 2000 K = (+ 


L) ? (K). 


(50) 


If the receiver aperture is given by a "top hat" function, 


I(K') becomes 


ERE? MEE 

J> G ae s: 

I(K') = K'D = MET ES RI 
2 2 





where b is given by 
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b ap ин 


Making these substitutions into (7.8) gives 
P_(K') = 21 3k? ($) 11/09 (K)I (K) Pia 


- f LdzC ^(z')cos(Kz) 


= OO 


Following Jarem, let 


A ÓN) GER 
giving 
Ели = (FE) è (K) I (K) [ azcy” (2) = [ dzcy”(2)cos (Kz) 
ZL (7.12) 


Il 

ryj 
X 

A 


(On Ex lower case n subscript indicates the dimensioned 
variable, while upper case N subscript indicates the 


dimensionless variable). Making similar substitutions into 


OB leads to 
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vhere 
В (ә) = = f azcy (2) / dKI (K) 6 (K) sin” (52) 
— c 0 
x Jg (aKl7 25) (7.14) 
and 
a = dipl . 


Another important relation is how the Hankel transform 
scales under the change of variables. Beginning with (7.4), 


nung (7.12) and (7.13) gives 
ә _ 122 2 
F (K) = 27 Ј в озу (ак o)d'edo . 
0 


cp 


P,(K) = SE | odoB,(o)Jyl(ak" 
0 


50. GIG) 


With the relations we have developed, (7.9) becomes 


2 


— 2 
FP. (K) = G (K) [2Cy = CN A BELG) 


where 


GR) = (K)I(K) . 


The solution to the integral equation becomes 


2 tier 2 iKz 
Cy Üz) = = E (K)e ZGK (Cais 
Or 
E -F. (K) 
2 el x = 2 IKZ 


This is the ill-posed solution of the integral equation. 
Tikhonov proposes using a low pass filter function q(K) such 
that the high spatial frequencies contaminated by the 
inversion process are removed. We should note that this 
process does nothing to errors at low spatial frequencies. 


The regularized solution is of the form 


— | иә | 
AZ) 2 J a" 5 E qikle aK. 719) 
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A profiling system would measure B, (0°). Using the 
relations developed, B (p') is scaled to В, (о) and trans- 


formed to F (K). In (7.19), all the quantities are knowns 


with the exception of e ES is given by 


N 
El 
Il 


2 
N İ aze, (2) ; 


which can be defined as a dimensionless form of Fried's ro 


by 
> 
Ro = 204^ è 
where 
22 4.6667 (Ҝ)-5/6 25/3 Gu. 
О sec 6 L 0 r i 


where $ is the zenith angle. Instruments exist to measure 
ro which immediately gives Ro: 

Another unknown, in terms of a practical measurement 
System, is L. In reality, the value of L need only to be 
Known approximately. The CHÉ profile can be normalized by 
independent measurements of 99 and rg. 

25”  .ıçting the Filter Function 
Clearly, in any measurement, noise will be present. 


Hence, the spectrum, FG, can be represented by a "true" 


component and a "noise" component, 
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N (K) 


7 , (7/21) 


F (K) = ES CHE, 


where the superscript t and n imply true and noise. Hence, 


Cy (EI is also composed of a "true" and "noise" component. 


If we let 
2 M 2 
Ci (K) = 2Cy Cu (K) ; 
(7.15) becomes 
F_(K) = G(K)Cz” (K) T 
" N : : 


The kernel, G(K), is type I to good approximation in 
Tikhonov's representation. Type I implies that there are 
no zeros along the real axis and as K + e, B(K) = О, 

For type I kernels, the filter function, q(K) can be 


represented by 


2 
"s ии (2.235 


GRO SE 2 
[G(K)]“ + M(K) 


where M(K) is a smoothing function. If the spectral 
characteristics of the noise are known, and the spectrum 


GA) unknown, M(K) becomes 


M(K) = axe? , (7.24) 
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where M(K) is a pan order stabilizer with regularization 


parameter a. If the noise spectra and signal CR 


spectra are known, one can construct an optimal Wiener 


filter with 





N(K) ' (7: 25) 


where S(K) is the signal spectra and N(K) is the noise 
spectra. 

If the noise process is white, like the shot noise 
in a measurement system, Tikhonov states that P = 1/2 in 
(7.24). Following Jarem, we can make an approximation to 
the optimal Wiener filter, (7.25). Consider an error 


functional of the form 
e = e oa - N | ale > (7.206) 


where o i2) has both signal and noise components and t in 


2 
the second term denotes the "true" component of CN (z). SBY 


Parseval's theorem, we can rewrite (7.26) as 


€ = 


> f EE GE 


ji. 





Now, GE is given by 


e E 2 i | ¡US "X 
C K = C _ — —— — — 5 
Ы ü G (K) [G(K)]^ * M(K) 


and ag is given by (7.16). Substituting (7.16), (7.28) 


and (7.21) into (7.27) gives 


Гә ик ии [G (K) 1 “ ` 
e aa J ¿Ey A — 2 IÓN 
Ea P [G(K)]^ * M(K) 


t 
F. (K) | 2 


E SENE dK ° (7.29) 


Putting (7.29) over a common denominator and performing the 


appropriate cancellations leaves 


` 00. 2 
x KW K , 
7: [16 (K)]% + m(x)]% 
(7.50) 
za 


vhere the second bracket in the numerator is merely CN (K). 


Squaring the numerator of (7.30) gives 


n 2 n ZE 
[G(K)F, (K)]* * 2M(K)G(K)F, (K)C4^ "(K) 


=>; | 


== С0 


[[(G(K)]* + M(K)]* 


[M (K) c, ^ (K) 1 ^] 
* — EDI RO aK (7556) 
11(6(K)17 + M(K)] 


L20 


darem notes that the middle term, 


n 2t 
2M (K)G (K) P. (K) Cy (K) 





Вер ии = 


since the noise spectra, P. (K) and the signal C ŻE 0%) are 


N 
produced by uncorrelated processes. Hence, the error 


becomes 


[IG (K) ,P (x)17 + [M (c, ^5 co 1? 


m 

|" 
N 
d 
— 8 


GR . (7.32) 
[[G(K)]* + M(K)]? 


We want to construct M(K) so that the error is minimized. 
This implies that the derivative of the error functional 


with respect to M(K) vanishes or 


de _ 
dM (K) = 0 a x 


After performing the differentiation, ve find that the 


minimizing M(K) is given by 
1 ee — : (7.34) 


: : n 
Since our noise process is white, BE (K) can be approximated 


: 2t K 
by a constant N. Jarem makes an assumption about Cy (K) 


lz 


that is invalid for vertical soundings. His data is along a 
2 km horizontal path. Hence, Cum Üz) is basically a 


constant making M) behave as 


A sin K 
K . 


Substituting back into (7.34) gives 


The average value of EE K is 1/2 making M(K) 


M(K) = ак“ È | (7.35) 
where 
a = ры : (6.26) 
A 
This corresponds to a stabilizer of order p = 1 and 


represents an optimal Wiener filter. It is interesting to 
note that M(K) is completely independent of the aperture. 
For vertical paths, Cu U (2) varies over five orders 
2t 


of magnitude. Hence, the assumption of a constant Cn (2) 


is not valid. However, it does represent a first order 
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approximation and is a starting point for numerical 
simulation. If problems occur in the inversion, one can 
model с̧ 2212) as an exponential decreasing in z, giving a 
new smoothing function M(K). Therefore, we have two 
Stabilizers of order p = 1 and p = 1/2 to attempt successful 
inversion. 

For Jarem's inversions, he finds that the regulari- 


E for the horizontal 


zation parameter, a, is nominally 10 
paths. Chapter VIII investigates this inversion technique 
on the 500 m vertical resolution data of Good shown in 
Figure 3.2. But before moving on to the simulations, I will 
investigate two analytical results that we will find 
convenient. 

In (7.14), the z integration can be completed 


analytically. First, interchanging the K and Z integrations 


glves 


2 K “z 
Ke 


В (о) = + [ aka Goa (ak1/29) | azcy?(2)sin 
0 — со 
x 


After completing the z integral, ve have 


co 


В (о) - i [f akxsqgoag(akl 75) (2657 - cy il, (7.38) 
0 
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where the sin” term is expanded using the double angle 
relation as before. For test purposes, it is convenient to 
5 : 2 

input cy” (2) as a cosine wave. Hence, Cy (K) becomes a 


delta function giving 


i İ aks (k)ə) (ak177,) tə6,? EX UE o 
0 


В, (ә) = 


2 


If the cosine wave has no offset, 25: İs zero and the K 


integral has an analytic result, 


Eu a = G(K')Jy(ak'"/*p) , (7.40) 


where K' is the wavenumber of the input cosine wave. Thus, 
(7.40) is an analytical check for softwave. 
3. Noise Computations 
Detectors in an instrumentation package that measure 
B, (ot) will add shot noise to the signal. From [Ref. 51], 


the average, rms shot noise current is given by 
... = 
iy (Wi Ze lao, (7:2 


vhere Av is the frequency bandvidth, I is the average 


current from the detector and e is the electron charge. I 
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depends on the intensity of the starlight and quantum 
2 clenoevaof the detectors. To find the irradiance of a 
given star on the telescope aperture, we can use the 


following [Ref. 52] relation, 


I 


1 
log M^! = . 4 (Magi - Mag5) , (7.42) 


where Ii, I, are the stellar intensities on earth and Magy 
and Mag» are the magnitudes of the different stars. Using 
the sun as a reference, we find that the star 9 Auriga 


8 


(Capella) has an irradiance of 1.05 x DNE W/m? on the 


telescope for Planck radiation between 0.4-0.6 um (visible). 


-10 
For the apodized aperture, the detector sees 2.15 x 10 W. 
The quantum efficiency is [Ref. 53] 
n a 
Q.E. - == , (7.43) 
q 


where n,- is the number of photoelectrons emitted by the 


detector/second and n, is the number of incident quanta/ 


q 
Second on the detector. Hence, Ng can be found using the 


following relation, 


ә 


(7.44) 


> 
Il 
-y 
< 
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where P is the power/wavelength interval of the starlight on 
the detector, h is Planck's constant and v is the frequency 
of the incident light. For the illustration of Capella, nq 
is 5.4 x 107 quanta/second given 2 = 500 nm. If we consider 
a detector with a quantum efficiency of 10%, l is 8.66 x 


1071” A. Hence, the variance of the current due to shot 


27 A” given a bandwidth of 500 kHz. 


noise is 1.39 x 10° 
Thus, we can construct a theoretical signal-to-noise ratio 
given by 
= 22107 
SNR 2 —l . -. $8meo x 10 . — 2 230 eni 


: 3.72 x 10 13 


We can use this figure for constructing Gaussian white noise 
in Chapter VIII. As pointed out earlier, the shot noise 


1/2 


increases with I Since the signal is proportional to 


the intensity, the signal-to-noise ratio increases with 
Ti E 
4. Inversion Conclusions 

The derivation of sections l and 2 is completely 
general in terms of the receiver apertures. Hence, the 
receiver apertures can be apodized to optimize inversion 
performance. We hope that the additional spatial filtering 
done by the aperture averaging will aid in regularizing the 


ill-posed nature of the inversion process. 
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VIII. NUMERICAL SIMULATION OF THE AMPLITUDE COVARIANCE 
FUNCTION INVERSION WITH TIKHONOV REGULARIZATION 
In this chapter, we hope to see if GER profiling by 

direct inversion of the scintillation amplitude covariance 
function with Tikhonov regularization is a practical 
measurement technique. We can accomplish this by 
numerically modeling the inversion process using the cig 
balloon data of Good presented previously. From this data, 
I construct the covariance function and use Tikhonov 
regularization to invert the covariance function regaining 
the original data. What follows is a step-by-step analysis 


of the numerical simulation. 


REC ONOTRUCIENG THE o: AND SCINTILLATION SPECTRA 

The finest step in obtaining gs is to transform Good's 
data into the dimensionless set of variables described in 
Chapter VII. We accomplish this by applying (7.11b) 
obtaining 570) from NE) cy” (K) is simply the Fourier 
transform of OR ZIE Por this test, the path length L is 
31.5 km representing the first 63 points of the known sei 
profile (500 m to 31.5 km). I obtained the value of Cy“ at 
the surface by assuming the -4/3 dependence of oe at low 


altitudes described in Chapter IV. Since the balloon data 


has a vertical resolution of 500 m, 4z = 0.0159. The spatial 
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frequency resolution of the Fast Fourier Transform (FFT) is 


given by [Ref. 54] as 


where N is the number of data points. Hence, AK is 0.9843. 
(We should remember that K and z are dimensionless 
variables.) The value of AK implies that NAT Should span 
zero to 63. However, FFT's are based on the Nyquist 
criteria, [Ref. 55] that requires two points per Az 
interval. Therefore, half the frequencies in the FFT are 
duplicates, implying that K spans zero to 31.5. 

The analysis in Chapter VII requires that Cine be an even 
function. This is convenient Since the Fourier transform of 
an even function is real. Figure 8.1 Shows the real part of 
C The modulus of RI is simply the absolute value 
of the real part since the imaginary part is zero. We 
obtain the scintillation spectrum F (K) by applying (7.16). 

It is interesting to note that the FFT is very robust in 
the presence of noise. Using an algorithm from [Ref. 56], I 
inJected Gaussian white noise with poorer signal-to-noise 
than that described in Chapter VIII, into the real and 
imaginary parts of CAUSE The real part of MONEO) with 
noise is shown in Figure 8.2. If we inverse FFT Ga" (EI with 


noise, we obtain the cy” (z) profile in Figure 8.3. The 
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noise has very little affect on the transform. We will find 
that very small amounts of noise produce drastic results in 


the inversion of the covariance equation. 


B. CONSTRUCTING THE COVARIANCE FUNCTION 

I determine the covariance function using (7.38). The 
numerical integration uses the Trapezoid rule at each 
detector separation p. We should note that the Trapezoid 
rule integrates functions with noise quite well since it 
does not require higher order derivatives in the data. 
Figure 8.4 shows the covariance function for a maximum 
detector separation of d = 0.5 m. The detector optics have 
a "top hat" aperture with a diameter of l cm. The form of 
the covariance function in Figure 8.4 is consistent with 
theoretical calculations [Refs. 6, 15]. 

The analytical form of the covariance function for a 
delta function cy” (K) is given by (7.40). 1 used this to 
check the covariance function software. The relative error 
is .02. Hence, the covariance function calculations have 
"noise" (errors) due to the numerical computations. We will 
use this as a source of noise for ill-conditioning. 

We should note that the covariance function has the 
proper dependence on optics aperture size and electro- 
magnetic wave wavelength. Figure 8.5 is the dependence of 
the covariance function on aperture size. As expected, the 


covariance function becomes wider with larger apertures. 


PSZ 


W 





UOTĄDUNĄJ QIUETIEAOI SpnąTTdwy UOTĄETTTĄJUTOS PSZTTEWAOŃ SUL р"8 әлпбтд 


G/MON 'NOIIBSBd3S 8012313ü G3ZI IBWeON 


Cc — 
5 

O m 
x 

D 

i 

— 
. IN 
Se m 
tJ 

O 
e 
С < 
D 

A 

tH 

> 
са 4 
m 


I 


W 2-30” S-UMHEMu”" 
10" =Y3J13W8BTIA SIIILHO ‘W G--NOLlLUMUd4S 40123130 ХЫЫ 
JONHIJHAOD Q3ZI IHWAON 


e 


9ZTS eanqazasdy sotado 1aAtaoday uo eouopuoedeq uorjounj SDUETIEAO) 


0/МОЈ "NOllUMUd38 4015231340 G(3ZI lBHWelON 





100-W S'= 'HSUH-W I'= 'QI10S-W 100"=WBIA SITIO 
W 1-382 S=HOBWH "NM S”-NOLllUMUd48 40123130 ХЫН 
ЗӘМӘТЈЫЛОЗӘ G3ZI IHWIAON 


G°8 sanbtTa 


JONBISBAOO ü3ZI IHWAON 


134 


This is predicted by the transform properties involved. 
Figure 8.6 shows that the covariance function widens for 
increasing wavelength. 


Sr profiling system based on this technique, we 


ın acC 
vould measure the intensity covariance function. Hence, the 
covariance function is the starting point for inversion to 
obtain the ə profile. The folloving sections describe the 
data manipulations required for an actual measurement 
system. It is interesting to note that experimental 
covariance function measurements of laser beam scintillation 
do not have the negative portion that represents the 
correlation of maxima and minima in the scintillation 
pattern [Ref. 28]. This may be due to the fact that the 
photomultipliers measure a positive definite quantity (i.e., 
intensity). 
C. CONSTRUCTING THE SCINTILLATION SPECTRUM AND APPLYING 

TIKHONOV REGULARIZATION 

The scintillation spectrum F (K) is the Hankel transform 
of the covariance function. The dimensionless form of the 
REGuStonmeis given by (7.15). 0 is found by (7.16). 
To compensate for the ill-conditioning present, ve introduce 
a lov pass filter of the form given in (7.23). This is the 
Tikhonov regularization. Figure 8.7 shovs CR) after the 


esch 3 
regularization with a stabilizer of the form ak . Note the 


affect of the filter on the higher spatial frequencies. 
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Figure 8.8 shows how EE changes with stabilizers of 
different powers of K. We obtain En (2) by inverse FFT. 
Figure 8.9 shows the true 202) profile, unregularized 


solution and regularized solution with a = 1 x l0 and a 


ҝ3 


stabilizer. Remember that a determines the amount of 
ill-conditioning allowed to enter the solution. The 
unregularized solution (dots) bears no resemblence to the 
original S profile. The regularized profile with a = 1 


MD a 


does very well resolving much of the 500 m 
resolution structure present in the actual data. The small 
"hump" in the regularized solution after the real data has 
gone to zero is part of the ill-conditioning. Note that the 
regularized solution has problems near zero. It in fact 
becomes negative (not shown) which is unphysical. This is 
very likely due to the frequency truncation. Consider a 
"top hat" aperture. We can Hankel transform the "top hat" 
to obtain the Airy spectrum.  Hankel transforming back to 
aperture coordinates produces an oscillation in the original 
aperture function, Figure 8.10. This is due to frequency 
truncation. This same truncation occurs in our calculation 
of the covariance function and a. Spectrum. Hence, an 
induced oscillation in cy” (K) transforms as the large 
negative value in the space coordinate, z. The fact that it 


is negative implies that the oscillation in the spectra is 


out of phase with the other frequencies present. 
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Lowering the regularization parameter a allows more of 
the ill-posed nature of the problem into the solution. 
Figure 8.11 shows the true solution and regularized solu- 


2 -13 3 stabilizer. 


tions with a = 1 x 10 ^ and 1 x 10 for agk 
Clearly, the choice of a is a trade-off between structure 
and magnitude accuracy. The aperture averaging does in fact 
filter the higher spatial frequencies, however manipulating 
a produces more favorable results than changing the optics 


diameter. 

It is interesting to note the affect on GRACZ) for 
stabilizers with different powers of K and the same 
regularization constant, a. Figure 8.12 shows cu (2) for 


K2, K? and K  stabilizers all Miente = o A K 


2 


Chapter VII, we found that K” stabilizers correspond to the 


: : : : 3 3 
optimal Wiener filter. From this simulation, I feel that K 


stabilizers offer the best fit to the c profile. 


D. REGULARIZATION CONCLUSIONS 


Initially, we were interested in whether direct 
inversion of the covariance function with Tikhonov 
regularization offered a practical way of obtaining cus 
profiles. I have shown that one can obtain high vertical 
resolution profiles (4z = 500 m) with this technique. 
However, covariance function measurements are very difficult 


to make. The quality of the measured covariance functions 


will never be comparable to those calculated numerically in 


142 


0 'siojeueadeq uorajezrde[nboyN 3uərzəjjytG Ҹдтм (2) 


30N111718 O3ZI IWWYON 
I 5” GE C SÌ S za Еј 


O xx 
Ee i 
bo 


100-£v4(21-30” 1)-544uM 'HSUBü-Ev»(6-30' 12-938 'üros-3neL 


ЗӘПЈӘЗЈӘ ИН 107 Н1Т1М (3193ANI 


(642 





PS INDIA 


Z 

O 

7 

x 

D 

[= 

— 

N 

M 

У 

Dcl o 

z 

PU 

N 

091 ~ 
022 

: QNO 


143 


N 


SST TABLES AUS TS TA OI 53 zI 8 Sm 


€ 


JAMATI JE UZZTIJEACJON 
C 3- S v Er 


















HSHT-bpvM()= 'IIJOS-EvMO- 'Ll0d-2vM(BHd 1B) S83ZI JIBHLS 


11-38 I=HHd1H '3efis3dBgB W IO' HLIM CILNIANI 


(2 INS 


° 
ә 


Ob 


08 


özl 


051 


auc 


ch 


(ZENO GAZI IBWelON 


144 


this section. Based on this fact alone, I feel that other 
techniques are more likely to yield accurate, high vertical 
resolution e ^ profiles. I will propose a method for 


obtaining such profiles in the next chapter. 
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IX. CONCLUSIONS 


The importance of the atmospheric isoplanatic angle and 
refractive turbulence profiles to the Space Defense 
Initiative (SDI) cannot be overstated. In this dissertation 
I have demonstrated an improved isoplanatic angle 
measurement system and presented a technique to measure 
refractive turbulence profiles. In the following sections, 


I will briefly review each area and propose future work. 


A.  ISOPLANATIC ANGLE MEASUREMENT 
We saw that the original Walters isoplanometer could be 
improved by apodizing the aperture of the receiver tele- 


— altitude veighting 


scope. To best approximate the z 
required for isoplanatic angle measurement, I constructed a 
double annulus aperture function optimized for a Celestron 
C8 telescope. The double annulus aperture weighting 
function is a significant improvement over that of the "top 
hat" clear aperture used on the previous isoplanometers. 
Sensitivity analysis of the double annulus weighting 
function to strong low altitude turbulence and inner scale 
changes with altitude showed no effect detrimental to a 
measurement system. Hence, we fabricated a fourth 


generation isoplanometer based on the double annulus 


aperture. This instrument is designed for night use and 
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offers many advantages over the previous system. First, the 
double annulus no longer requires the use of the Celestron 
Cl4 telescope as did the "top hat" aperture. This is a 
considerable savings in cost and instrument managability. 
This instrument also has the ability to make isoplanatic 
angle measurements once per second giving data with very 
high temporal resolution. Trends in the isoplanatic angle 
are emerging now that have never been seen before and 
indicate that this instrument may be an extremely effective 
tool to probe the tropopause. 

In Chapter V, I detailed tests proving that the 
instrument is indeed measuring the isoplanatic angle 
correctly. And lastly, I presented data from Maui showing 
interesting trends in the isoplanatic angle. Currently, two 
fourth generation isoplanometers are operating at different 
locations, supporting high priority programs. 

Clearly, the greatest need in the area of isoplanatic 
angle measurement is a painstaking correlation between the 
structure observed in isoplanatic angle measurements and 
meteorological data. One hopes that the isoplanometer might 
effectively measure turbulence at the tropopause. This 
would save money in balloon flights, etc. 

Also of interest would be to modify the fourth 
generation system to provide high temporal resolution data 


during the day. This requires the use of high quality 
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telescope mounts and drives, but I feel daytime data is well 


worth the expense. 


B. REFRACTIVE TURBULENCE PROFILING 

In this section, we investigated the inversion of the 
scintillation amplitude covariance equation. I included the 
filtering done by the receiver optics and obtained a closed 
form solution to the inversion problem. Unfortunately, 
inverting the covariance function is ill-posed in a 
mathematical sense, meaning that noise in the data causes an 
unstable solution. Hence, I incorporated a regularization 
technique developed by Tikhonov to compensate for the ill- 
conditioning. Tikhonov regularization is a low pass filter 
designed to eliminate high spatial frequencies most 
contaminated by the ill-conditioning. 

I modeled the inversion numerically using an actual 
refractive turbulence profile with 500 m vertical 
resolution. We saw that the Tikhonov regularization 
reproduced much of the profiles very well. However, we note 
that actual covariance function measurements are very 
difficult to make. Hence, I recommend pursuing the 
weighting function approach to refractive turbulence 
profiling described in Chapter VII. One should be able to 
obtain suitable weighting functions given an appropriate 
basis set. This would allow refractive turbulence profiling 
by scintillation intensity variance measurements, operating 
much like the isoplanometer. 
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APPENDIX A 


OTHER WEIGHTING FUNCTIONS 


Many different weighting functions can be generated 
using the technique mentioned in Chapter III. What follows 
is a brief overview of the many weighting functions I 


generated while trying to approximate to the 2213 


behavior 
or find weighting functions that peak at a given z. 

First, a few words about the double annulus aperture. 
In Figure 3.8, we have noted the error remaining at low z. 
This can be improved by making the inner annulus larger, 
Figure A.1. However, you sacrifice some performance at 


larger altitudes. Because of the 27/3 


weighting for the 
isoplanometer, we felt error should be minimized at higher z 
at the expense of low z. 

I also tried other combinations of Bessel functions. 
One is the "wedding cake" aperture function (in my 
nomenclature), which is composed of two "top hat" functions 
with different transmissivities, Figure A.2. Typically, the 
"wedding cake" aperture gave a weighting function of z 
below 3 km, Figure A.3. Associated with the "wedding cake" 
is the "three tier" aperture function, which has one more 


"layer". Weighting function behavior with the "three tier" 


is basically the same as the "wedding cake". 
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One can incoherently add spectra of different aperture 
functions since the correlation distance is so small (on the 
order of centimeters), given a telescope with a sufficiently 
large aperture. Using this technique, I was able to make 
the weighting function bend at a given z, Figure A.4. The 
weighting function in Figure A.4 is an incoherent addition 
of the "top hat" and "wedding cake" aperture functions. 
Weighting functions with this characteristic are also 
important. 

From the Bateman manuscripts, I tried some Hankel 
transform pairs. If we consider the Hankel transform of 


f(x), (the spectrum), in the form 


со 


Е(у) = İ (300 (xy) (xy) 1/2ax (1) 
0 
the function х71/22-ах has the spectrum — T zyj RKA 


The exponential aperture function has the weighting function 


shown in A.5. This weighting function is also on the z? 


side of 27/3, 

Shown in Figures A.6 and A.7 are the weighting functions 
associated with sine and cosine spectra. Both have 
analytical aperture functions that can be found in the 
Bateman Manuscript [Ref. 34]. 

This is by no means a complete list of all the apertures 
I found weighting functions for. However, these are the 


more interesting and applicable. 
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APPENDIX B 
MAUI METEOROLOGICAL BALLOON DATA 
LO-CATE WO-8000 UPPER AIR DATA MANAGEMENT SYSTEM P106277 


REVZ.11 (C) 1984 
DATE(GMT) TIME(GMT) ASCENT NO. FLIGHT EQ. STATION SONDE NO. 


850917 04:16 17 5 0 6685 
TIME PRESSURE HEIGHT TEMP DP-DEP RH 
MIN. MBARS M-MSL CELA CEL. PERCENT 

0.00 711.4 0 12.0 zp 2256 

0.65 700.0 135 Yə 17.6 28.3 

Ig 75 682. 4 348 KEE 14.6 67 

2.48 659. 4 633 9.5 15.0 34.7 

4. 00 GET 1038 4.6 1280 41.7 

6.10 514,0 1709 0.0 yarar 42.9** 

6.15 576.4 Jo; -0.1 11.3 42.6 

6.28 572 1788 90.8 17.8 46.9 

6.35 517 1816 -0.3 EEN 43.3 

7. 80 540.6 2261 -4.1 2 sən 

575 518.8 2585 -5.5 12.6 36.8 

8.82 sak 2608 -2.7 1226 3:5 

8.88 51577 2632 -5.9 1 37.6 

9.67 500.0 2874 -7.3 1351 34.6 

9.88 495.7 2941 -7.9 1550 29.0 
11.40 466.0 3418 -11.6 HB 38.4 
12.80 443.3 3800 -12.6 14.5 28.9 
1356 400.0 4576 -18.4 14.0 28.3 
—” 305.6 6530 -32.5 10, 5 55 
275 300.0 6660 -33.5 1050 35.6 
25 260.8 7627 -41.6 6.6 47.9 


158 


TIME 

MIN. 

26.60 
767 
29.80 
31330 
34, 78 
356.72 
5713 
57.40 
38.40 
20.45 
75:56 
40.38 
43.53 
44.83 
45.30 
45.63 
48.10 
52.60 
54.15 
56.53 
20:75 
59.42 
00705 
56.17 
60.73 
66.75 
68.03 
70.00 
70.07 
70. I3 
70. 25 


PRESSURE 
MBARS 


250.0 
233:8 
2161 
200.0 
IO DNE 
15222 
1 500 
148.0 
140.6 
140.1 
T3970 
126.9 
109.0 
UNC 
100.0 
98.4 
87.4 
70.0 
64.9 
5 
51.6 
50.0 
41.6 
3060 
20.0 
7: 
10.0 
0.0 
0.0 
0.0 
030 


HEIGHT 
M-MSL 


7913 

8359 

8876 

5777 
10427 
11022 
11147 
11227 
11531 
11554 
11606 
12143 
T3029 
3396 
13527 
13620 
14310 
15633 
16092 
16830 
17495 
17690 
18835 
20883 
23427 
23554 
27711 
28794 
28794 
28794 
28794 


LSD 


TEMP 


SELER 


“2536 
-47.5 
“20.2 
254.2 
25520 
E) 
140/596 
z / 1.31 
—””: 
-26.4 
“0:1 
“2:7 
“7: 
"76.6 
12057 
bs 
— 72.2 
— c 
- 65. 0 
-65.0 
-62.0 
OIOM 
=59. 4 
Say 
EDO 
EDU 
5005 
“69.3 
-- 
-67.6 
757: 


DP-DEP 


CEE: 


RH 


PERCENT 


38.4 
5757 
5055 
7:5: 
o 
34.8 
34.7 
34.7 
34. 7 
39.4 
34.7 
34.4 
34.3 
34.2 
34.1 
34.1 
34.5 
34.8 
34.9 
34.9 
34.9 
34.9 
34.7 
34.4 
34.4 
34.4 
34.2 
34.2 
34.2 
34.3 
34.4 


TIME 
MIN. 
702 
71.78 
71785 
/ 1.92 
71.98 
12403 
72.10 


TIME 
MIN. 


On DM A UU N El O 


il 
12 
13 
14 
15 
16 
17 
18 
19 


PRESSURE 
MBARS 


HEIGHT 
M-AS 


204 

424 

718 
1038 
1358 
1677 
2015 
2329 
2668 
2978 
3292 
3582 
3856 
4139 
4421 
4708 
5000 
5292 
55095 


HEIGHT 


M-MSL 


28794 
30992 
30992 
30992 
30992 
30992 
30992 


TEMP 
CECS 
-64.7 
-43.0 
-42.0 
-41.0 
-40.3 
=3 974 
—55 


WIND VELOCITY 


DIRECTION 
DEGREES 


79 
102 

99 

94 

49 

43 

67 
116 
178 
188 
216 
247 
246 
212 
T99 
198 
206 
219 
234 
220 
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DP-DEP 


CEL. 


SPEED 
KNOTS 


2020 
LTS 
11.6 
6.8 
25 
3.4 
Ste 
333 
555 
e 
Gare 
Gan 
ә 
9.4 
12. 
10.7 
9. 
10.4 
rz. 
12. 


RH 
PERCENT 

34.4 
391 
38.4 
39.0 
59.7 
40.7 
41.8 


HETGHT 
FT-MSL 
0.0 
668.7 
139m. 1 
2556975 
3405.5 
44538 
550254 
6612925 
7641.8 
8753.4 
9769.3 
1 080m2 
117500 
12652283 
135781 
14503.9 
15444.7 
16404.0 
173622 
1832272 


TTME HEIGHT DIRECTION SPEED 


MIN. M-AS DEGREES KNOTS 
20 5577 234 10.6 
ZI 5159 255 Jo 1 
22 6462 2506 105 
7 5770 248 217: 
24 7083 246 207 
25 57 265 21 1 
26 HIS 284 2550 
27 8054 292 3303 
28 5395 295 575: 
29 8662 296 77 
30 8943 ZE 55275 
sal 9277 283 51557 
32 9588 268 25, 6 
33 9889 259 25002 
34 101901 240 2752 
55 10494 2 28.0 
36 10801 229 3727 
57 ll 07 246 30.0 
38 11409 260 29.0 
7 J.L 7 30 284 20.1 
40 5022 247 23.4 
41 0 243 2 
42 nes 248 18.4 
43 12879 235 Pe” 
44 Fri 251 J 0 
45 13443 209 IS 
46 1 5725 197 14.4 
47 14002 203 050 
48 14282 24 2:8 
49 14575 50 LORE 
50 14868 19 ДӘ. 


17: 


RETERE 
FT-MSL 
192816 
20240.38 
2120070 
2221025 
23239.4 
24268.4 
27116 
26423.3 
2754159 
28418.8 
2000 
30433738 
5155057 
32445.6 
51:70 
34428.0 
201077 
36440.2 
374532. 
38485.2 
39464.0 
40408.2 
412210 
4225006 
430782 
44103.5 
4502157 
45939. 2 
40555 
47816.8 
45751. 4 


TIME 
MIN. 
SI 
52 
53 
54 
55 
56 
5: 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 


HEIGHT 
M-AS 


15163 
15457 
15751 
16048 
163535 
16665 
16970 
17270 
17568 
17874 
18189 
18504 
ığı” 
19459 
20116 
2077 
24364 
27603 
28794 


DIRECTION 
DEGREES 


87 
90 
102 
107 
95 
80 
86 
97 
100 
DM 
88 
89 
96 
106 
106 
106 
106 
106 
27 
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SPEED 
KNOTS 


19% 
15m 
20.4 
18.6 
1920 
20.4 
2319 
2328 
Za 
1070 
Saa 
H972 
25.4 
22 
27298 
27m 
275 
211 
47.8 


HEIGHT 
FESHMSE 
49746.0 
50 7190885 
5167738 
52649.4 
5365657 
54674.6 
5507527 
56660.0 
5 763/89. 
58641.1 
5969759 
6076/81009 
61742.8 
63842.6 
6599825 
68154. 3 
7055 
90560. 9 
94468.3 
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